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^ ' Abstract 



We present a study of D = 4 supersyminetric Yang-Mills matrix models with S0{3) 



^ . mass terms based on the cohomological approach and the Monte Carlo method, 

^s^ , In the bosonic models we show the existence of an exotic first/second order transition 

On I from a phase with a well defined background geometry (the fuzzy sphere) to a phase with 

commuting matrices with no geometry in the sense of Connes. At the transition point the 
sphere expands abruptly to infinite size then it evaporates as we increase the temperature 



m 

o 

CN ' (the gauge coupling constant) . The transition looks first order due to the discontinuity in 

the action whereas it looks second order due to the divergent peak in the specific heat. 

The fuzzy sphere is stable for the super symmetric models in the sense that the bosonic 

p\ , phase transition is turned into a very slow crossover transition. The transition point is 

j^ , found to scale to zero with N. We conjecture that the transition from the background 

sphere to the phase of commuting matrices is associated with spontaneous supersymmetry 

breaking. 

The eigenvalues distribution of any of the bosonic matrices in the matrix phase is 
found to be given by a non-polynomial law obtained from the fact that the joint proba- 
bility distribution of the four matrices is uniform inside a solid ball with radius R. The 
eigenvalues of the gauge field on the background geometry are also found to be distributed 
according to this non-polynomial law. 

We also discuss the D = 3 models and by using cohomological deformation, localization 
techniques and the saddle-point method we give a derivation of the D = 3 eigenvalues 
distribution starting from a particular D = 4 model. 



*An early version of this work titled " Noncomniutative Gauge Theory As/Is Matrix Models Around Fuzzy 
Vacua And Emergent Geometry" was presented to BM Annaba University as a habilitation thesis on 20 January 
2011. 
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1 Introduction 

It can be argued using the principles of quantum mechanics and classical general 
relativity that the picture of spacetime at the very large as being a smooth manifold must 
necessarily break down at the Planck scale Xp [HIE]- At this scale localization looses its 
operational meaning due to intense gravitational fields and formation of black holes and as 
a consequence one expects spacetime uncertainty relations which in turn strongly suggest 
that spacetime has a quantum structure expressed by [j;^,Xj^] = iXpQ^u- The geometry 
of spacetime at the very small is therefore noncommutative. 

Noncommutative geometry [23], see also [73H78] and [79], allows for the description of 
the geometry of smooth differentiable manifolds in terms of the underlying C*— algebra of 
functions defined on these manifolds. Indeed given the following three data a) the algebra 
A = C°°{M) of complex valued smooth and continous functions on a manifold M, b) the 
Hilbert space Ti = L?'{M, S) of square integrable sections of the irreducible spinor bundle 
over M and c) the Dirac operator V = ^^{d^j^+^oj^ahl"'!^) associated with the Levi-Civita 
connection uj one can reconstruct completely the differential geometry of the manifold 
M. These three data compose the so-called spectral triple {A^T-L^V) corresponding to 
the Riemannian manifold M . In the absence of spin structure it is sufficient to use the 
Laplacian A instead of the Dirac operator in the spectral triple. 

Noncommutative geometry is also more general than ordinary differential geometry in 
that it also enables us to describe algebraically the geometry of arbitrary spaces (which 
a priori do not need to consist of points) in terms of spectral triples. The paradigm 
of noncommutative geometry adopted so often in physics is to generalize the ordinary 
commutative space M by replacing the commutative algebra ^ by a noncommutative 
algebra Aq. The result of this deformation is in general a noncommutative space Mq 
defined precisely by the spectral triple {Ag^T-ie^Vg/ /S.0) where the Hilbert space T-Lq is the 
representation space of the noncommutative algebra Aq and Vq/ iS.Q is the deformation of 
the commutative Dirac operator/Laplacian P/A |23y79j. 

Noncommutative geometry was also proposed (in fact earlier than renormalization) 
as a possible way to eliminate ultraviolet divergences in quantum field theories [8OII8T]. 
The quantum spacetime of [ll[2] is Lorentz-covariant based on the commutation relations 
[x^, Xy] = iX^Q^u with Q^^ satisfying [xx, Q,,u] = 0, Q^.uQ'"' = and {le^^xpQ'^'Q^^f = 
1. As it turns out quantum field theory on this space is ultraviolet finite [82j which 
is a remarkable consequence of spacetime quantization. This result is however not very 
surprising. Indeed this phenomena of "regularization by quantization" already happens in 
quantum mechanics. For example while classical mechanics fail to explain the blackbody 
radiation in the ultraviolet (the UV catastroph) quantum mechanics reproduces the correct 
(finite) answer given by the famous experimentally verified Stefan-Boltzman law. 

Noncommutative field theory is by definition a field theory based on a noncommutative 
spacetime |26U27j . The most studied examples in the literature are the Moyal-Weyl spaces 
Rg which correspond to the case Q/^i^ = ^^i/ where 9^,^ are rank 2 (or 1) antisymmetric 
constant tensors, i.e. 

[Xfj,,Xu] = i9ij,u- (1-1) 



This clearly breakes Lorentz symmetry. The corresponding quantum field theories are 
not UV finite |125j and furthermore they are plagued with the so-called UV-IR mixing 
phenomena |83j . 

For example in the case of scalar field theory on Moyal-Weyl spaces the UV-IR mixing 
was shown to destroy the perturbative renormalizability of the theory in |lU9yilU] . The 
UV-IR mixing in this model means in particular that the non-planar contribution to the 
two-point function which is finite for generic values of the external momentum behaves 
as l/iOpY' in the limit p — > and/or 9 — > 0. The physics at very large distances is 
thus altered by the noncommutativity which is supposed to be relevant only at very short 
distances. Equivalently it was shown in jlll] that although the renormalization group 
equations are finite in the IR regime their perturbative approximations are IR divergent. 
In [112H114] a modified scalar field theory on noncommutative R^ was proposed and shown 
to be renormalizable to all orders in perturbation theory in d = 2 and 4. The modification 
consists of making the action invariant under the duality transformation of [115] by adding 
a harmonic oscillator potential to the kinetic term which precisely modifies the free theory 
as desired. 

For scalar field theory with (j)'^ interaction the phase diagram consists of three phases 
instead of the usual two phases found in commutative scalar field theory. We observe a 
disordered phase, a uniform ordered phase and a novel nonuniform ordered phase which 
meet at a triple point possibly a Lifshitz point |116pil7] . The nonuniform phase is the 
analogue of the matrix phase in pure gauge models (see below) in the sense that in this 
phase the spacetime metric is modified by quantum fiuctuations of the noncommutative 
field theory since the Laplacian at the transition point is found to be (9^)^ and not 
dfi |116) . In the nonuniform ordered phase we have spontaneous breakdown of translational 
invariance |117H119] . The transition from the disordered phase to the nonuniform ordered 
phase is thought to be first order and the nonuniform ordered phase is a periodically 
modulated phase which for small values of the coupling constant is dominated by stripes 
|117] . This should hold even in two dimensions. 

This behavior was confirmed in Monte Carlo simulations on the noncommutative torus 
[10H12] in [Unilm] and on the fuzzy sphere [SUED] in p^2til2i] . 

It is therefore natural to conjecture that there must exist two fixed points in this 
theory, the usual Wilson-Fisher fixed point at = and a novel fixed point at ^ = oo 
which is intimately related to the underlying matrix model structure of the model |126j . 
The crucial input in a scalar field theory on a noncommutative space is the scalar potential 
which when written in the matrix base is a random matrix theory which for a (j)^ interaction 
is given by 

V = N{aTr]Vf + hTrW^). (1.2) 

The main property of this matrix model which plays a central role in the phase structure 
of a noncommutative (f)'^ interaction is the existence of one-cut/two-cut transition, with 
the one-cut (disordered) phase corresponding to a>a^: where |132pi33j 

a, = -2A (1.3) 



The kinetic term is trying to add a geometry to the dynamics of the matrix M which is 
at the heart of the rich phase structure we observe. For small b the usual Ising model 
transition is expected and the = fixed point should control the physics. We expect 
on the other hand that the 6 = oo fixed point should control most of the phase diagram 
since generically the cp'^ interaction, i.e. the coupling b is not weak. 

Noncommutative gauge theories attracted a lot of interest in recent years because of 
their appearance in string theory [12 1 184 ^ 185]. For example it was discovered that the 
dynamics of open strings moving in a fiat space in the presence of a non- vanishing Neveu- 
Schwarz B-field and with Dp-branes is equivalent to leading order in the string tension 
to a gauge theory on a Moyal-Weyl space Rg. Extension of this result to curved spaces 
is also possible at least in one particular instance, namely the case of open strings mov- 
ing in a curved space with S^ metric. The resulting effective gauge theory lives on a 
noncommutative fuzzy sphere S^ [86fl88j . 

Another class of noncommutative spaces, besides Moyal-Weyl spaces, which will be 
very important for us in this article is fuzzy spaces |89y9Uj. Fuzzy physics is by definition 
a field theory based on fuzzy spaces [9T][92] . The original idea of discretization by quanti- 
zation (fuzzification) works well for co-adjoint orbits such as projective spaces. A seminal 
example of fuzzy spaces is the fuzzy two-dimensional sphere S^ |29ll30j . The fuzzy sphere 
is defined by three N x N matrices Xi, i = 1,2,3 playing the role of coordinates operators 
satisfying ^^ xf = 1 and the commutation relations 

1 N"^ -1 
[xi,Xj] = iOeijkXk , 9 = —— , C2 = . (1.4) 



The fuzzy sphere which is the simplest among fuzzy projective spaces was actually pro- 
posed as a nonperturbative regularization of ordinary quantum field theory in |127j . See 
also |98lll01ill021ll28] . However we know now that this can not be correct because of 
the complicated phase structure of noncommutative (/>^ scalar field theories on the fuzzy 
sphere [12911130] . Instead it is becoming very clear that it is more appropriate to think 
of the fuzzy sphere as a nonperturbative regularization of noncommutative quantum field 
theory on Moyal-Weyl plane. For example we can show in a particular double scaling 
limit that the nonperturbative phase structure of scalar fields on the Moyal-Weyl plane 
can be rigorously identified with the nonperturbative phase structure of scalar fields on 
the fuzzy sphere, and in particular the matrix or nonuniform ordered phase on the fuzzy 
sphere goes precisely to the stripe phase on the Moyal-Weyl plane [131] . 

In 3— dimensions we have the fuzzy three sphere S^ [93]. In 4— dimensions we have 3 
fuzzy manifolds which are obtained from co-adjoint orbits. The direct product of fuzzy 
two spheres S^ x S^ [94] , fuzzy GV\ [95|l96] and fuzzy "8,% as a squashed fuzzy CP^ [97] . 

The most appealing aspect of discretization by quantization remains its remarkable 
success in preserving continuous symmetries including supersymmetry and capturing cor- 
rectly topological properties j3ll[9ll[92] . Indeed the fuzzy approach does not suffer from 
fermion doubling [9814100] . it extends naturally to supersymmetry |101H104] and it cap- 
tures correctly nontrivial field configurations such as monopoles and instantons using the 
language of projective modules and K-theory [105H108] . 



The noncomniutative Moyal-Weyl spaces should be thought of as infinite dimensional 
matrix algebras not as continuum manifolds. Quantum fluctuations of a gauge theory on 
Rg will generically make the vacuum which is here the Moyal-Weyl space itself unsta- 
ble. To see this effect we formally rewrite gauge theory on R^ as a matrix model with 
A^— dimensional matrices D^ where A^ — > oo. This is given by 

252 



S= ^-::^' Trn[^[D.,b,]--Br^] . (1.5) 



By computing the effective potential in the configuration Di = —(f)B^- Xj we will verify 
explicitly that the Moyal-Weyl space itself, i.e the algebra [x^,Xu] = iO^^, ceases to exist 
above a certain value g* of the gauge coupling constant. We find g'^ = tt'^/N in four 
dimensions (see section 2). In other words by moving in the phase diagram from strong 
coupling to weak coupling the geometry of the Moyal-Weyl spaces (including the star 
product and the representation of operators by fields) emerges at the critical point (7* . As a 
consequence a natural regularization of Moyal-Weyl spaces must include in a fundamental 
way matrix degrees of freedom. The regularization which will be employed in this article is 
given by fuzzy projective spaces [134j . The main reason behind this choice is the fact that 
the phenomena of emergent geometry |20pi35pi36j which we observed here on Moyal-Weyl 
spaces shows up also in all matrix models on fuzzy projective spaces. 

Let us note the deep connection which exists between the geometry in transition we 
observe here on Moyal-Weyl space R^ and the perturbative UV-IR mixing and beta func- 
tion of U{1) gauge theory computed in |137lll38j . The structure of the effective potential 
in 2 dimensions indicates that there is no transition and no UV-IR mixing in U{1) gauge 
theory on R^ and therefore the theory is expected to be renormalizable. This was shown 
numerically to be true in [139j . We also expect that noncommutative M = 1 supersym- 
metric U{1) gauge theory in 4 dimensions is renormalizable. 

Let us note also that fuzzy regularization is different from the usual one which is based 
on the Eguchi-Kawai model [13] and the noncommutative torus |40fl42j . i.e the twisted 
Eguchi-Kawai model. The twisted Eguchi-Kawai model was employed as a nonperturba- 
tive regularization of noncommutative gauge theory in [140H142] where the instability and 
the phase transition discussed here were also obtained. 

Finally we note that the strong coupling phase (above g^,) of noncommutative gauge 
theory on Rg corresponds to 9 = 00. It is dominated by commuting operators. The 
limit 9 — > 00 is the planar theory (only planar graphs survive) |125j and it is intimately 
related to large A^ limits of hermitian matrix models [14311144] and [117| . In this phase of 
commuting operators supersymmetry may be broken. 

The phenomena of emergent geometry associated with noncommutative gauge theory 
is therefore a major motivation behind choosing fuzzy projective spaces as a nonpertur- 
bative regularization of Moyal-Weyl spaces. The other motivation comes of course from 
the novel phase known variously as stripe, nonuniform ordered or matrix phase found in 
noncommutative scalar field theory which we have only briefly discussed in this introduc- 
tion. The phenomena of emergent geometry is also, on the other hand, intimately tied to 
reduced Yang-Mills models which will be the topic of central interest in the remainder of 
this introduction and in most of this article. 
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It is well established that reduced Yang-Mills theories play a central role in the non- 
perturbative definitions of M-theory and superstrings. The BFSS conjecture [3] relates 
discrete light-cone quantization (DLCQ) of M— theory to the theory of N coincident DO 
branes which at low energy (small velocities and/or string coupling) is the reduction to 
+ 1 dimension of the 10 dimensional U{N) supersymmetric Yang- Mills gauge theory [3]. 
The BFSS model is therefore a Yang-Mills quantum mechanics which is supposed to be the 
UV completion of 11 dimensional supergravity. As it turns out the BFSS action is nothing 
else but the regularization of the supermembrane action in the light cone gauge [5] . 

The BMN model [6] is a generalization of the BFSS model to curved backgrounds. It is 
obtained by adding to the BFSS action a one-parameter mass deformation corresponding 
to the maximally supersymmetric pp-wave background of 11 dimensional supergravity. 
See for example [THS]. We note, in passing, that all maximally supersymmetric pp-wave 
geometries can arise as Penrose limits of AdSp x S'^ spaces [TQ| . 

The IKKT model [Hj is, on the other hand, a Yang-Mills matrix model obtained by 
dimensionally reducing 10 dimensional U{N) supersymmetric Yang-Mills gauge theory to 
+ dimensions. The IKKT model is postulated to provide a constructive definition of 
type II B superstring theory and for this reason it is also called type IIB matrix model. 
Supersymmetric analogue of the IKKT model also exists in dimensions d = 3, 4 and 6 
while the partition functions converge only in dimensions d = 4,6 [52y53j. 

The IKKT Yang-Mills matrix models can be thought of as continuum Eguchi-Kawai 
reduced models as opposed to the usual lattice Eguchi-Kawai reduced model formulated 
in [13]. We point out here the similarity between the conjecture that the lattice Eguchi- 
Kawai reduced model allows us to recover the full gauge theory in the large N theory and 
the conjecture that the IKKT matrix model allows us to recover type II B superstring. 

The relation between the BFSS Yang-Mills quantum mechanics and the IKKT Yang- 
Mills matrix model is discussed at length in the seminal paper [12] where it is also shown 
that toroidal compactification of the D-instanton action (the bosonic part of the IKKT 
action) yields, in a very natural way, a noncommutative Yang-Mills theory on a dual 
noncommutative torus [24J. From the other hand, we can easily check that the ground 
state of the D-instanton action is given by commuting matrices which can be diagonalized 
simultaneously with the eigenvalues giving the coordinates of the D-branes. Thus at tree- 
level an ordinary spacetime emerges from the bosonic truncation of the IKKT action while 
higher order quantum corrections will define a noncommutative spacetime. 

In summary, Yang- Mills matrix models which provide a constructive definition of string 
theories will naturally lead to emergent geometry |20j and non-commutative gauge theory 
|145[|146] . Furthermore, non-commutative geometry [23p76j and their non-commutative 
field theories [26p27j play an essential role in the non-perturbative dynamics of superstrings 
and M-theory. Thus the connections between non-commutative field theories, emergent 
geometry and matrix models from one side and string theory from the other side run deep. 

It seems therefore natural that Yang-Mills matrix models provide a non-perturbative 
framework for emergent spacetime geometry and non-commutative gauge theories. Since 
non-commutativity is the only extension which preserves maximal supersymmetry, we 
also hope that Yang-Mills matrix models will provide a regularization which preserves 



sup er symmetry 

In this article we will explore in particular the possibility of using IKKT Yang-Mills 
matrix models in dimensions 4 and 3 to provide a non-perturbative definition of emergent 
spacetime geometry, non-commutative gauge theory and supersymmetry in two dimen- 
sions. From our perspective in this article, the phase of commuting matrices has no 
geometry in the sense of Connes and thus we need to modify the models so that a geom- 
etry with a well defined spectral triple can also emerge alongside the phase of commuting 
matrices. 

There are two solutions to this problem. The first solution is given by adding mass 
deformations which preserve supersymmetry to the flat IKKT Yang-Mills matrix models 
|14j or alternatively by an Eguchi-Kawai reduction of the mass deformed BFSS Yang-Mills 
quantum mechanics constructed in |15 p64f[66] . The second solution, which we have also 
considered in this article, is given by deforming the flat Yang-Mills matrix model in Z? = 4 
using the powerful formalism of cohomological Yang-Mills theory [32H34^ [68] . 

These mass deformed or cohomologically deformed IKKT Yang-Mills matrix models 
are the analogue of the BMN model and they typically include a Myers term [TBj and 
thus they will sustain the geometry of the fuzzy sphere [29l[30] as a ground state which at 
large N will approach the geometry of the ordinary sphere, the ordinary plane or the non- 
commutative plane depending on the scaling limit. Thus a non-perturbative formulation 
of non-commutative gauge theory in two dimensions can be captured rigorously within 
these models [I9l[l471[l48] . See also [TOlEl]. 

This can in principle be generalized to other fuzzy spaces ^31j and higher dimensional 
non-commutative gauge theories by considering appropriate mass deformations of the flat 
IKKT Yang- Mills matrix models. 

The problem or virtue of this construction, depending on the perspective, is that in 
these Yang- Mills matrix models the geometry of the fuzzy sphere collapses under quantum 
fluctuations into the phase of commuting matrices. Equivalently, it is seen that the geom- 
etry of the fuzzy sphere emerges from the dynamics of a random matrix theory |56p58j. 
Supersymmetry is naturally expected to stabilize the spacetime geometry, and in fact 
the non stability of the non-supersymmetric vacuum should have come as no surprise to 
us [35] . 

We should mention here the approach of [38] in which a noncommutative Yang-Mills 
gauge theory on the fuzzy sphere emerges also from the dynamics of a random matrix 
theory. The fuzzy sphere is stable in the sense that the transition to commuting matrices 
is pushed towards infinite gauge coupling at large A'' [57] . This was achieved by consider- 
ing a very special non-supersymmetric mass deformation which is quartic in the bosonic 
matrices. This construction was extended to a noncommutative gauge theory on the fuzzy 
sphere based on co-adjoint orbits in [39]. 

Let us also note here that the instability and the phase transition discussed here were 
also observed on the non-commutative torus in [72ll731[T^0Ul42j where the twisted Eguchi- 
Kawai model was employed as a non-perturbative regularization of non-commutative 
Yang-Mihs gauge theory [40ti42] . 

In this article we will study, using cohomological matrix theory and the Monte Carlo 



method, the mass deformed Yang-Mills matrix model in D = 4 as well as a particular 
truncation to D = 3. We will derive and study a one-parameter cohomological deformation 
of the Yang- Mills matrix model which coincides with the mass deformed model in Z) = 4 
when the parameter is tuned appropriately. We will show that the first/second order 
phase transition from the fuzzy sphere to the phase of commuting matrices observed in 
the bosonic models is converted in the super symmetric models into a very slow crossover 
transition with an arbitrary small transition point in the large A^ limit. We will determine 
the eigenvalues distributions for both D = 4 and D = 3 throughout the phase diagram. 
The D = 3 eigenvalues distribution can be obtained from a particular Z? = 4 model by 
means of the methods of cohomological Yang-Mills matrix theory, large N saddle point 
and localization techniques [361137]. 

This article is organized as follows: In section 2 we give a brief discussion of the 
phenomena of emergent geometry on the Moyal-Weyl space using the effective potential. 
In section 3 we give a short review on fuzzy projective spaces. In section 4 we review 
results on emergent geometry in bosonic D = 3 Yang-Mills matrix models. In section 5 
we will derive the mass deformed Yang-Mills quantum mechanics from the requirement 
of supersymmetry and then reduce it further to obtain Yang- Mills matrix model in D = 
4 dimensions. In section 6 we will derive a one-parameter family of cohomologically 
deformed models and then show that the mass deformed model constructed in section 5 
can be obtained for a particular value of the parameter. In section 7 we report our first 
Monte Carlo results for the model D = 4 including the eigenvalues distributions and also 
comment on the D = 3 model obtained by simply setting the fourth matrix to 0. In 
section 8 we solve a particular D = A model and show by means of cohomological Yang- 
Mills matrix theory and localization techniques that it is equivalent to the matrix D = 3 
Chern-Simons theory. By using the saddle point method we derive the D = 3 eigenvalues 
distribution. We conclude in section 9 with a comprehensive summary of the results and 
discuss future directions. The detail of the simulations are found in appendices A and B 
while appendix C contains a detail calculation of the super symmetric mass deformation 
and appendix D contains a calculation of the star product on the Moyal-Weyl plane 
starting from the star product on the fuzzy sphere. 

2 Gauge Theory on Noncommutative Moyal-Weyl 
Spaces 

The basic noncommutative gauge theory action of interest to us in this article can be 
obtained from a matrix model of the form (see [26] and references therein) 



V^'^det(vri?)^^ ^2 _ VO'detiTTB)^^ /^^ ^ ^ 1 



S = ^'TT'-' Tr^F- = y"-2-' Trn ^[D., D,] - ^B;^ • (2.1 



2g' "^ '' 2g 



«j 



2 



Here i,j = 1, ..., d with d even and 9 has dimension of length squared so that the connection 
operators Di have dimension of (length)"^. The coupling constant g is of dimension 
(mass) ~ 2 and B^^ is an invertible tensor which in 2 dimensions is given by B^- = e~- = 



-Cij while in higher dimensions is given by 



B- 



( 



\ 



tjj 



(2.2) 



-*i / 



The operators Ai belong to an algebra A. The trace is taken over some infinite dimensional 
Hilbert space l-L and hence rr-^[l)j, Z)j] is /O in general, i.e. Try\Di^Dj\ is in fact a 
topological term |12j . Furthermore we will assume Euclidean signature throughout. 
Minima of the model (j2.ip are connection operators Di = Bi satisfying 



m,Bj 



^B.\ 



(2.3) 



We view the algebra ^ as ^ = Mat„(C) (E> An- The trace Try^ takes the form Try^ = 
TrnTr-y^ where "Hn is the Hilbert space associated with the elements of ^„. The config- 
urations Di = Bi which solve equation ()2.3p can be written as 



Bi 



]B^^^ 



(2.4) 



The operators Xi which are elements of An can be identified with the coordinate operators 
on the noncommutative Moyal-Weyl space R^ with the usual commutation relation 



[Xii Xj\ 



iOB, 



ij- 



Derivations on Rg are defined by 



Indeed we compute 



di = iBi. 



[di,Xj] 



H]- 



(2.5) 



(2.6) 



(2.7) 



The sector of this matrix theory which corresponds to a noncommutative U{n) gauge field 
on Rg is therefore obtained by expanding Di around Bi(^ln- We write the configurations 



Di 



--B^.^Xj®ln + Ai, i+ 



Ai. 



(2.^ 



The operators Ai are identified with the components of the dynamical U{n) noncommu- 
tative gauge field. The corresponding U{n) gauge transformations which leave the action 
(|'2.ip invariant are implemented by unitary operators U = exp(iA) , UU^ = U^U = 



1 , A+ = A which act on the Hilbert space Ti = Tin 



Tin as Di — >UDiU~^, i.e. 



Ai — >UAiU^ — iU[di, U~^] and Fij — >UFijU~^. In other words U{n) in this setting must 



be identified with U{%n 



S 



) %n)- The action ()2.ip can be put into the form 

^Je<^det{TTB) 



W 



-Trnsn 



C\2 
ij ) ■ 



(2.9) 
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The curvature F[^ where C is a U{n) index which runs from 1 to n^ is given by 



^3 

1 

2 



F§ = [4, if] - [d„Af] - tf^Bc{Af.Af} + :-dABc[At,Af]. (2.10) 



In calculating F^- we used [Ta,Tb\ = ifABcTc, {Ta,Tb} = (IabcTc and TtTaTb = ^. 
More explicitely we have defined Ta = -f- for the SU{n) part and Tq = ^=1„ for the 



C/(l) part. The symbols cIabc ^^^ defined such that date are the usual SU{n) symmetric 

symbols while dabo = daob = doab = y^^^ab, daoo = and dooo = y f- 

Finally it is not difficult to show using the Weyl map, viz the map between operators 
and fields, that the matrix action (|2.9|) is precisely the usual noncommutative U{n) gauge 
action on R^ with a star product * defined by the parameter 9Bij [2U1[57]- In particular 
the trace Tr-}^^ on the Hilbert space Tin can be shown to be equal to the integral over 
spacetime. We get 



^g'J 



d^x [F^^f , i^f = a,Af - d.A't - ^fABc{Af,Af}, + '-dABc[At,Afl. 

(2.11) 



Let us note that although the dimensions dim?^ and dimT^^ of the Hilbert spaces Ti and 
Tin are infinite the ratio dim'H/dim7^„ is finite equal n. The number of independent 
unitary transformations which leave the configuration (j2.4p invariant is equal to dim?^ — 
dim?^^ — n?. This is clearly less than dimH for any n > 2. In other words from entropy 
counting the U{1) gauge group (i.e. n = 1) is more stable than all higher gauge groups. 
As we will show the U{1) gauge group is in fact energetically favorable in most of the 
finite A^ matrix models which are proposed as non-perturbative regularizations of (|2.ip in 
this article. Stabilizing U{n) gauge groups requires adding potential terms to the action. 
In the rest of this section we will thus consider only the U{1) case for simplicity. 

Quantization of the matrix model (|2.1|) consists usually in quantizing the model (|2.11|) . 
As we will argue shortly this makes sense only for small values of the coupling constant 
g^ which are less than a critical value g^. Above g'^ the configuration Bi given by ()2.4p 
ceases to exist, i.e. it ceases to be the true minimum of the theory and as a consequence 
the expansion (j2.8p does not make sense. 

In order to compute this transition we use the one-loop effective action which can be 
easily obtained in the Feynamn-'t Hooft background field gauge. We find the result 

T = S + ^TuTr^d In (v^6ij - 2iT^j^ - Tr^^ \nV\ (2.12) 

The operators V^ = 'D{Di , Pj and Fij act by commutators, viz 2?^(..) = [l)j, [Z)j, ..]], 
T)i[..) = [Di, ..] and J-ij{..) = [Fij,..]. Next we compute the effective potential in the 
configuration Z)j = —(f)B^- Xj. The curvature Fy in this configuration is given by 9 Fij = 
{O'^cfP' — i)B^- . The trace over the Hilbert space Ti is regularized such that Tr-^1 = A^ is 
a very large but finite natural number. We will also need J2i j ^ij ^ij = d. The effective 
potential for d ^ 2 is given by 

a{9'^(t>^ -1)'^ + In (j). (2.13) 



(d - 2)A2 
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The coupling constant a is given by 

- 
"=— 2- X^ ■' = ''-'"■ P'") 

We take the hmit N — > oo keeping X^N fixed. It is not difficult to show that the 
minimum of the above potential is then given by 



- ^i^. (2.15) 

The critical values are therefore given by 

a. = l^A^iV = ^^. (2.16) 

Thus the configuration Z)j = —(j)B~- Xj exists only for values of the coupling constant 
A which are less than A*. Above A* true minima of the model are given by commuting 
operators, i.e. 

i[Bi,Bj]=0. (2.17) 

By comparing with p.3p we see that this phase corresponds to = oo. The limit 9 — > oo 
is the planar theory (only planar graphs survive) |125] which is intimately related to large 
N limits of hermitian matrix models |117j . 

This transition from the noncommutative Moyal-Weyl space ()2.3p to the commuting 
operators (|2.17p is believed to be intimately related to the perturbative UV-IR mixing [83] . 
Indeed this is true in two dimensions using our formalism here. 

In two dimensions we can see that the logarithmic correction to the potential is ab- 
sent and as a consequence the transition to commuting operators will be absent. The 
perturbative UV-IR mixing is, on the other hand, absent in two dimensions. Indeed in 
two dimensions the first nonzero correction to the classical action S in the effective action 
(I2.12P is given by 






J- — b I ^ad ,T->9 •' ij ^T->2 *i + ••• 



Jk '^ JpV 

S + 2jTr\F,,ip)\'J^^j-^{l-cose,,p,k,)\. (2.18) 



ip 

,2 



By including a small mass m and using Feynman parameters the planar and non-planar 
contributions are given respectively by 



k"^ + m? {p — k)"^ + m? 4tt Jq z^ 
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In above z is defined by z^ = (OijPi)'^ {m? + x(l — x)p^) and Ki(z) is the modified Bessel 
function given by 



oo 2 ^2 



Z ZG 

zKi{z)= I dte-*e-tF = l+ In -- + .... (2.21) 

Jo 2 2 

We observe that in two dimensions both the planar and non-planar functions are UV finite, 
i.e. renormalization of the vacuum polarization is not required. The infrared divergence 
seen when m? — > cancel in the difference 11^ — 11^^. Furthermore 11^ — II'^^ vanishes 
identically in the limit 9 — > or p — > 0. In other words there is no UV-IR mixing in 
the vacuum polarization in two dimensions. 

The situation in four dimensions is more involved |137pi38] . Explicitly we find that the 
planar contribution to the vacuum polarization is UV divergent as in the commutative 
theory, i.e. it is logarithmically divergent and thus it requires a renormalization. It is 
found that the UV divergences in the 2—, 3— and 4— point functions at one- loop can be 
subtracted by a single counter term and hence the theory is renormalizable at this order. 
The beta function of the theory at one-loop is identical to the beta function of the ordinary 
pure SU{2) gauge theory. The non-planar contribution to the vacuum polarization at one- 
loop is UV finite because of the noncommutativity and only it becomes singular in the 
limit of vanishing noncommutativity and/or vanishing external momentum. This also 
means that the renormalized vacuum polarization diverges in the infrared limit p — >0 
and/or 9 — > which is the definition of the UV-IR mixing. 

We expect that supersymmetry will make the Moyal-Weyl geometry and as a conse- 
quence the noncommutative gauge theory more stable. In order to see this effect let A^, 
a = 1, ..., M be M massless Majorana fermions in the adjoint representation of the gauge 
group U{T-L). We consider the modification of the action (|2.ip given by 



S' = S+:^^^^^Y.TrnXa,ADuXa]. (2.22) 



The irreducible representation of the Clifford algebra in d dimensions is s = 22 dimen- 
sional. Let us remark that in the limit 9 — > the modified action S has the same 
limit as the original action S. By integrating over \a in the path integral we obtain the 
Pfaffian (pf(7jPj)) . We will assume that pf(7jPj) = (det(7jPj))^. The modification of 
the effective action (j2.12p is given by 



TrsTr^A In (v^ - ^-lHj^iA . (2.23) 



It is not very difficult to check that the coefficient of the logarithmic term in the effective 
potential is positive definite for all M such that Ms < 2(1 — 4. For Ms = 2d — A the 
logarithmic term vanishes identically and thus the background (j2.4p is completely stable 
at one-loop order. In this case the noncommutative gauge theory (i.e. the star product 
representation) makes sense at least at one-loop order for all values of the gauge coupling 
constant g. The case Ms = 2d — Am. d = A (i.e. M = 1) corresponds to noncommutative 
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N" = 1 supersymmetric ^(1) gauge theory. In this case the effective action is given by 

r' = S + -TrdTr^A In ( 5ij - 2i—^Fij j - —TrgTr^A In f 1 - -lHj^J^ij \ ■ (2.24) 

This is manifestly gauge invariant. In 4 dimensions we use the identity Tvg'jijjjk^i = 
s{SijSki — 6ikdji + dii6jk) and the first nonzero correction to the classical action S is given 
by the equation 

,d-2 



S+{ g I)7'''adp2-^«J2?2-^«J 



5 + 2(1 



d-2. 



Tr\Fi,{p)\' 



k A;2 {p - k) 



;{l- cos eijPikj)\. (2.25) 



This correction is the only one-loop contribution which contains a quadratic term in the 
gauge field. The planar and non-planar corrections to the vacuum polarization are given 
in this case by 



n^ 



1 



1 



1 



fcA;2(p-/c)2 (47r)i 70 



dx 



dt 



t2 



2-1 



~x{l—x)p^t 



(2.26) 



n 



NP 



1 1 

kk'^ {p- k)' 



■ cos OijPikj 



\ f^ dt 
dx / — ; — e 



(47r)2 Jo Jo t2 ^ 



■x(l-x)p2i_teL 



(2.27) 



The planar correction is UV divergent coming from the limit t — > 0. Indeed we compute 
(including also an arbitrary mass scale /x and defining 6 = 2—2) 



n^ 



-^fdx{p-)i-\x{l-x)t,f^--T{2-^) 
(47r)2 Jo ^^ 2 



{^^ 



2\-e 



(47r)2 



1 fl p2 

7— / dxlnx(l — x)— ^ + 0(e) 

e Jo /^^ 



(2.28) 



The singular high energy behaviour is thus logarithmically divergent. The planar correc- 
tion needs therefore a renormalization. We add the counter term 



SS = -2(1 - l^)i^l / d'^xF? 



(4^) 



u 



-2(l-^)^-/l4(p)|^. (2.29) 
o (47r)2 



The effective action at one-loop is obtained by adding ()2.25p and the counter term (12.290 . 
We get 

1 



252 (^) 



\Mp)Y 



(2.30) 



252 (^) 



1 
1 



2+2(1 



d-2 



m' 



n 



NPn 



2(1 



d-2 (/i2)-^l 



(4vr)^ 



+ 



3 1 



2^2 2(47r) 



■7 



dx\nx{l 






-n 



NP 
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(2.31) 



This equation means that the gauge couphng constant runs with the renormahzation scale. 
The beta function is non-zero given by 

/3(5(/^))=/^^ = -Y|^/(/i). (2.32) 

The non-planar correction is UV finite. Indeed we compute the closed expression 

n'"'' = TTxi [ dxKoiz) , z^ = {e,jp^fx{l - x)p\ (2.33) 

In the limit 9 — > and/or p — > we can use Kq[z) = — In | and obtain the IR singular 
behaviour 



1 r...J%p^)Mi-xy 

(47r)2 



n^^^ = _^^/ dxln ''^"^'' ' "^ . (2.34) 



In summary although the Moyal-Weyl geometry is made stable at one-loop order by the 
introduction of supersymmetry we still have a UV-IR mixing in the quantum gauge theory. 
The picture that supersymmetry will generally stabilize the geometry will be confirmed 
nonperturbatively in this article whereas the precise connection to the UV-IR mixing 
remains unclear. 



3 Fuzzy Spaces 



3.1 The Fuzzy Sphere 

The ordinary sphere S"^ is defined in global coordinates by the equation 

3 



nGR3,J]n2 = l. (3.1) 

a=l 

A general function can be expanded in terms of spherical harmonics 

oo 
f{n) = ^fkmXkmin). (3.2) 

fc=0 

Global derivations are given by the generators of rotations in the adjoint representation 
of the SU{2) group, namely 

^a = —ieabcnbdc , [Ca,Cb] = iCabc^c- (3.3) 

The Laplacian is given by 

C = LaCa , eigenvalues /(/ + !) , / = 0, ..., oo. (3.4) 

According to [79 tll97j all the geometry of the sphere is encoded in the K-cycle or spectral 
triple [A.U.C?)- -A = C°°(S'2) is the algebra of ah functions / of the form ([SJ]) and 
T-L is the infinite dimensional Hilbert space of square integrable functions on which the 
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functions are represented. In order to encode the geometry of the sphere m the presence of 
spin structure we use instead the K-cycle {A, Ti, T>, 7) [23] where 7 and T> are the chirahty 
and the Dirac operators on the sphere. Remark in particular that 

Yiioc rii +in2 , liooc ns jYi-ioc ni - in2. (3.5) 

The fuzzy sphere is a particular deformation of the above triple which is based on the fact 
that the sphere is the co-adjoint orbit SU{2)/U{1), namely 

9(^39^^ = riaCTa , g£SU{2) , n^S"^. (3.6) 

This also means that S"^ is a symplectic manifold and hence it can be quantized in a 
canonical fashion by simply quantizing the volume form 

CO = sinOdOAdcj) = -eabcnaduhAdric- (3.7) 

The result of this quantization is to replace the algebra C°^ (S"^) by the algebra of matrices 
MatL+i. 

Mati+i is the algebra of (-^+1) x (L+1) matrices which acts on an (L+1)— dimensional 
Hilbert space Hl with inner product {f,9) = ■j^Tr{f^g) where f,g^MatL^i. The spin 
s = L/2 IRR of SU(2) is found to act naturally on this Hilbert space. The generators are 

[La, Lf,] = ieabcLc , J] ^a = C2 = ^(^ + !)• (3-8) 

a 

Generally the spherical harmonics Ifcm(^) become the canonical SU{2) polarization ten- 
sors Ykm [198J. They are defined by 



[La, [La,Yim]] = l{l + l)Yi^ , [L±,YiJ^ = ^/ {l^m){l±m + l)Yi^±i , [L3,Yi^] = mYi^(3.9) 
They also satisfy 

Yil = i-irYi.m , j^TrYi,^,Yi,ra, = {-ir'5i,i,5m,,-m,. (3.10) 

Matrix coordinates on S^ are defined by the k = 1 tensors as in the continuum, namely 

O Q Q - - t I J fj , •, 

X]^ + X2 + Xo = 1 , \Xa, XfiJ = — -^Gabc^c, Xa = — ;^. (3-11) 

VC2 y/C2 

"Fuzzy" functions on S\ are linear operators in the matrix algebra while derivations are 
inner defined by the generators of the adjoint action of SU{2), i.e. 



Cai<P) = [La,(P]. (3.12) 

A natural choice of the Laplacian operator A on the fuzzy sphere is therefore given by 
the following Casimir operator 

AL = Cl = [La,[La,..]]. (3.13) 
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Thus the algebra of matrices Mat^+i decomposes under the action of the group SU{2) as 

-®- = 001020. .eL. (3.14) 

The first L/2 stands for the left action of the group while the other L/2 stands for the 
right action. It is not difficult to convince ourselves that this Laplacain has a cut-off 
spectrum of the form /(/ + 1) where Z = 0, 1, ..., L. As a consequence a general function on 
S^ can be expanded in terms of polarization tensors as follows 

L I 
1=0 m=-l 

The fact that the summation over / involves only angular momenta which are <L originates 
of course from the fact that the spectrum l{l + 1) of the Laplacian A is cut-off at I = L. 
The commutative continuum limit is given by L — )-cxd. This is the first limit of interest 
to us in this article. Therefore the fuzzy sphere is a sequence of the following triples 

{MatL+i,HL,AL). (3.16) 

Again in the presence of spin structure the fuzzy sphere S|^ will be defined instead by the 
K-cycle (A, H, D, T) where the chirality operator F and the Dirac operator D are given 
for example in |100j and references therein. 

3.2 Fuzzy CP^ 

In this section we will give the K-cycle associated with the classical Kahler manifold 
CP (which is also the co-adjoint orbit SU{3)/U{2)) as a limit of the K-cycle which 
defines fuzzy (or quantized) CP when the noncommutativity parameter goes to 0. We 
will follow the construction of |134j . 

Let Ta, a = 1, ...,8 be the generators of SU{3) in the symmetric irreducible represen- 
tation (n, 0) of dimension N = 2(71 + l){n + 2). They satisfy 

[Ta,Tb] = ifabcTc . (3.17) 

Ta = ^nin + 3) = \n\^ , dabcTaTb = ^^^T,. (3.18) 

3 D 

Let ta = Aa/2 (where Xa are the usual Gell-Mann matrices) be the generators of SU{3) 
in the fundamental representation (1,0) of dimension A^ = 3. They satisfy 

'^tah = -Sab + {dabc + ifabc)tc 

trstah = -Sab , tr^tatbtc = jidabc + ifabc)- (3.19) 

The A— dimensonal generator Ta can be obtained by taking the symmetric product of n 
copies of the fundamental 3— dimensional generator ta, viz 

Ta = (ia0l0---0l + l(E>ta0...(E>l + ... + l0l«)...0ia)symmetric- (3.20) 
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In the continuum CP^ is the space of ah unit vectors |^ > in C^ modulo the phase. Thus 
e \ip >, for aU 0G[O, 2it[ define the same point on CP . It is obvious that all these vectors 
e IV' > correspond to the same projector P = |^ >< V'l- Hence CP is the space of all 
projection operators of rank one on C^. Let Hjv and H3 be the Hilbert spaces of the 
SU{?>) representations (n, 0) and (1,0) respectively. We will define fuzzy CP^ through 
the canonical SU{3) coherent states as follows. Let n be a vector in R^, then we define 
the projector 

P3 = h + nata (3.21) 

The requirement P3 = P3 leads to the condition that n is a point on CP satisfying the 
equations 

2 4 2 

[na,nb\ = , na = - , dabcnarib = -^c- (3.22) 

We can write 

P3 = |n,3>< 3,n|. (3.23) 

We think of |n, 3 > as the coherent state in H3 (level 3x3 matrices) which is localized at 
the point n of CP . Therefore the coherent state \n,N > in H^v (level N x N matrices) 
which is localized around the point n of CP is defined by the projector 

P/V = \n, N >< N, n\ = [Ps^Pa^-^Ps) symmetric- (3.24) 

We compute that 

1 Tl 

trstaPs =< n,3\ta\n,3 >= -Ua , tr^TaPN =< fi,N\Ta\n,N >= -Ua- (3.25) 

Hence it is natural to identify fuzzy CP^ at level N = |(n + l)(n + 2) (or CP^) by the 
coordinates operators 

Xa = -Ta. (3.26) 

n 
They satisfy 



2i 4 3 2 3 

[Xa.Xb] = —fabcXc , X^ = -(1 + -) , dabcXaXb = ^(1 + T^)^c- (3.27) 

n 6 n 6 In 

Therefore in the large N limit we can see that the algebra of Xa reduces to the continuum 
algebra of n^. Hence Xa — >na in the continuum limit A^ — t-oo. 

The algebra of functions on fuzzy CP^ is identified with the algebra oi NxN matrices 
MatN generated by all polynomials in the coordinates operators Xa- Recall that N = 
2{n + l)(n + 2). The left action of SU{3) on this algebra is generated by (n, 0) whereas 
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the right action is generated by (0, n). Thus the algebra Matj\f decomposes under the 
action of SU{3) as 

(n,O)0(O,n) = ®^=o(P,P)- (3-28) 

A general function on fuzzy CP^ is therefore written as 



p=0 



i{P:P) 



The Tp J y are SU{3) matrix polarization tensors in the irreducible representation {p,p)- 
I^,Is and Y are the square of the isospin, the third component of the isospin and the 
hypercharge quantum numbers which characterize SU{3) representations. 

The derivations on fuzzy CP^ are defined by the commutators [T^, ..]. The Laplacian 
is then obviously given by An = [Ta, [Ta, ■■■]]■ Fuzzy CP^ is completely determined by 
the spectral triple CP^ = {Mat^, A^v, Hat). Now we can compute 



n 

tTNFPN =< n, N\F\n, N >= F^(n) = J^ ^Xk^'^/VsV^'')- ^^'^^^ 

p=0 



Ap,p) 



The Yp J yin) are the SU{3) polarization tensors defined by 

Y^^!'jly{n) =< n, iVlT^ty In, iV > . (3.31) 

Furthermore we can compute 

trN[Ta,F]PN =< n,N\[Ta,F]\n,N>= {CaFN){n) , Ca = -ifabcnbdc- (3.32) 

And 

tTNFGPN =< n, N\FG\n, N >= Fn * G7v(n). (3.33) 

The star product on fuzzy CPjy is found to be given by (see below) 

Fn * GN{n) = ^ ^-Y^Ka^b^...Kapb^da^...da^FN{n)db^...dh^GN{n) 

2 

Kab = l^Sab - riaUb + {clabc + ifabcjric- (3.34) 



3.3 star Products on Fuzzy CP' 

The sphere is the complex projective space CP . The quantization of higher CP ~ 
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with A^ > 2 is very similar to the quantization of CP , and yields fuzzy CP r . In this 



section we will explain this result by constructing the coherent states |190H192] . the Weyl 
map [193] and star products p4HT96] on all fuzzy CP^~^ following pi] . 
We start with classical CP defined by the projectorqU 



1 n /2(A^-1) 

P = -l + ajvnX, «JV = -Y ^ ^ ' . (3.35) 

The requirement P^ = P will lead to the defining equations of CP as embedded in 
R ~^ given by 

2 1 ^ 2 N-2 
f^a = 1 ' dabcnanb = ^^"-c- (3-36) 

The fundamental representation N of SU{N) is generated by the Lie algebra of Gell-Mann 
matrices ta = -f,CL = 1,..., N"^ — 1. These matrices satisfy 

'^tah = JiT^ab^ + {dabc + "ifabcjic 

Trtahtc = -{dabc + ifabc) , Trtah = -^ , Trta = 0. (3.37) 

Let us specialize the projector (|3.35p to the "north" pole of CP given by the point 
no = (0,0, ..., 1). We have then the projector Pq = jjl + aNtj\[2_i = diag(0, 0, ..., 1). So at 
the "north" pole P projects down onto the state l-^o >= (0,0, ..., 1) of the Hilbert space 
C^ on which the defining representation of SU{N) is acting . 

A general point n£ CP can be obtained from uq by the action of an element 

g£SU{N) as n = gfiQ. P will then project down onto the state \Tp >= glipo > of C^ . 
One can show that 

P = 1^ X ^1 = 51^0 X ^o|5+ = gPQg+. (3.38) 

gtN2_^g+ = n^ta. (3.39) 

This last equation is the usual definition of CP^^^ . Under g — >gh where h^U{N — 1) 
we have htj^2_ih'^ = tjv^-i j i-e U{N — 1) is the stability group of tj^2_i and hence 

Q-pN-i ^ SU{N)/U{N - 1). (3.40) 

Thus points n of CP^~"^ are then equivalent classes [g] = [gh], h^U{N — 1) . 

The fuzzy sphere S^ is the algebra of operators acting on the Hilbert space Hi which 
is the (2s + 1)— dimensional irreducible representation of SU{2). This representation can 
be obtained from the symmetric product of L = 2s fundamental representations 2 of 
SU{2). Indeed given any element g^SU{2) its s— representation matrix U^^'{g) can be 
obtained as follows 

C/W(c,) = C/(2)(5)0^...0^[/(2)((,)^ 2s -times. (3.41) 



^In this section we use N to denote the dimension of the space instead of the size of the matrix approximation. 
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U^^\g) is the spin 1/2 representation oi g^SU{2) . 

Similarly fuzzy CP^ ~ is the algebra of operators acting on the Hilbert space Hi 
which is the Oh -dimensional irreducible reprsentation of SU{N). This dimension is given 
explicitly by 

d(^) - ( (^I^1±M (3 42) 

"'^ ~^{N-iy.{2s)\- ^^^ ^ 

Also this ds -dimensional irreducible representation of SU{N) can be obtained from the 
symmetric product of L = 2s fundamental representations N of SU{N). 

Remark that for s = 1/2 we have d,,J = N and therefore H^J = C^ is the fundamen- 
tal representation of SU{N). Clearly the states iV'o > and \ip > of H^^ will correspond 

in Hs to the two states \no,s > and \n,s > respectively so that iV'o >= Ko) ^ > and 
\il) >= |n, ^ >. Furthermore the equation |^ >= gli^o > becomes 

\n,s>=U^'\g)\no,s> . (3.43) 

U^^'{g) is the representation given by 

C/W(5) = U^^\g)(g)s...(^sU^'^\g),2s- times. (3.44) 

To any operator F on Hs (which can be thought of as a function on fuzzy CP^ ~ ) we 
associate a "classical" function Fs{n) on a classical CP by 

Fs{n) =< n, s\F\n, s > . (3.45) 

The product of two such operators F and G is mapped to the star product of the corre- 
sponding two functions 

Fs*Gs{n)=<n,s\FG\n,s> . (3.46) 

Now we compute this star product in a closed form. First we will use the result that any 
operator F on the Hilbert space Hs admits the expansion 

F= [ dix{h)F{h)U^'\h). (3.47) 

JsU{N) 

U^^'{h) are assumed to satisfy the normalization 

TrU^'\h)U^'Hh') = df^6{h'^ - h'). (3.48) 

Using the above two equations one can derive the value of the coefficient F{h) to be 

Fih) = ^TrFU^^Hh-'). (3.49) 

d)s 

Using the expansion (|3.47p in (|3.45|) we get 

Fsin)= [ dfi{h)F{h)J'\n,h) , J'\n,h) = < n,s\U^'\h)\n,s >. (3.50) 

JSU{N) 
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On the other hand using the expansion (j3.47p in (j3.46p will give 



F, * GJn) 



dfi{h)dfi{h')F{h)G{h')J'\n,hh'). 



(3.51) 



SU{N) 



The computation of this star product boils down to the computation of uj^'{n, hh ). We 
have 



J'\n,h) = <n,s\U^'\h)\n,s> 



^ 1, 

< n,-\(g)s... 



< n, ■ 



t/(^)(/i)®,...®,?7(N)(/i) 



1 
2 

(3.52) 



J^\n,h) = < ij\U^^\h)\iP > . 



(3.53) 



In the fundamental representation N of SU{N) we have U^ '{h) = exp(ir7T,"ta) = c(m)l + 
is"'{m)ta and therefore 

oj^2\n,h) = <i)\c{m)l + is''{m)ta\'ip >= c{m) +is''{rn) <'4)\ta\'ip > . (3.54) 



a;(i)(n,/i/j') = < ^1^^^^^^')!^' > 

= < V'|(c(m)l + is"-{ni)ta){c{rn)l + is" {ni )ta)\i) > 

= c{m)c{'m ) + i[c{m)s"'{m ) + c{m )s'^{m)] < Tp\ta\ip > 

- s"-{m)s^{m') < i^ltatbl'ip > ■ 

Now it is not difficult to check that 

<V'|t,|V'> = TrtaP=^n'' 

<tp\tatb\tp> = TrtahP = —5ah + ^{dahc + ifahc)n''. 



Hence we obtain 



u}^2'{^n,n) = c[m) + I s{ra).n. 



(3.55) 



(3.56) 



(3.57) 



L0^'2'[n,hh) = c[m)c[m ) —s[m).s[m)+i 



cimjs [m 



ON 



+ c{m )s"'{m) 
These two last equations can be combined to get the result 



n" - -fidabc + ifabc)n's''{ni)s\m ). (3.58) 



cj(2)(n,/t/t')_^(2)(n,/i)t^(2)(n,/i') = -—g(^rn).s{m)- -^{dabc + ifabc)n''s°'{m)s''{m') 



+ 



a 



N ^a^b 



rf'n Sa{m)sb{m ). 



(3.59) 
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We can remark that in this last equation we have got ride of ah reference to c's. We would 
like also to get ride of all reference to s's. This can be achieved by using the formula 

Saim) = —-J-2\n, h). (3.60) 

lUN on"- 

We get then 

J^^\n,hfi)-uj^^\n,h)uj^^\n,li) = Kab-?—OJ^^\n,h)-—^uj^^\n,}i). (3.61) 

2 1 

^a.b = -ZT^^ab - naUb -\ {dabc + ifabc)n''. (3.62) 



Fs*G,in) = Y,-^^-—^Ka,b^::Ka,b, 



Therefore we obtain 

2 

kl{2s-k)\ 

SU{N) on an fc 

(3.63) 

We have also the formula 

(3.64) 
This allows us to obtain the final result |134j 

2s 



k=0 



I.— n "^ ' 



Specialization of this result to the sphere is obvious. In the last appendix we will discuss 
a (seemingly) different star product on the fuzzy sphere which admits a straightforward 
flattening limit to the star product on the Moyal-Weyl plane. 

Let us also do some examples. Derivations on CP^~^ are generated by the vector 
fields Ca = —ifabcnbd/dnc which satisfy [Ca,Cb\ = ifabc^c- The corresponding action on 
the Hilbert space Hg is generated by La and is given by 

< n,s\U^'\h-^)FU^'\h)\n,s>=< no,s\U^'\g-'^h-'^)FU'^'\hg)\no,s > . (3.66) 

U^^'{h) is given by U^^\h) = exp{ir]aLa). Now if we take r] to be small then one computes 

< n,s\U^''\h)\n,s >= 1 + irja < n,s\La\n,s > . (3.67) 
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On the other hand we know that the representation U^^'{h) is obtained by taking the 
symmetric product of 2s fundamental representations N of SU{N) and hence 

< n, s\U^'\h)\n, s >= (< n, ^|1 + ir/,t,|n, ^ >f' = 1 + ^(2s)7?„^n,. (3.68) 

In above we have used the facts La = ta(E>s----^sta, \n,s >= |n, ^ > (^jg-.-^sln,^ > and 
the first equation of (j3.56p . Finally we get the important result 

< n, s\La\fi, s >= sa^na- (3.69) 

We define the fuzzy derivative [La,F] by 

{CaF)s{n) = <n,s\[La,F]\n,s> 



(3.70) 



= SON 


Ua * Fs{n) - Fs *na{n) 


= ifabcrf Q^^Fsin). 


Finally we note the identity 


-l-TrFG= [ Fs*Gs{n). 



(3.71) 



4 Review of Bosonic D = 3 Yang-Mills Matrix 
Models 

The principal goal of the present article is the construction of a new nonperturbative 
method for noncommutative gauge theories with and without supersymmetry based on 
a class of Yang-Mills matrix models in which the classical minima are given by fuzzy 
projective spaces. These matrix models are directly related to the celebrated IKKT matrix 
model |li yi45[ri46yi49| in d = 3 and d = 4 dimensions with mass deformation which may 
or may not preserve supersymmetry. 

The flat IKKT models, i.e. without mass deformation, are obtained by dimensionally 
reducing U(A^) super Yang-Mills theory in flat d dimensions onto a point, i.e. to zero 
dimension. The dynamical variables are d matrices of size N with action 

S = -^Tr[X^,X,]^ + Tri^T^[X^,i;]. (4.1) 



The partition functions of these models are convergent in dimensions d = 4, 6, 10 |48H54j . 
In d = 3 the partition function may be made finite by adding appropriate mass deformation 
consisting of a positive quadratic term in the matrices X^ which damps flat directions. In 
d = 3, 4 the determinant of the Dirac operator is positive deflnite |46y50j and thus there 
is no sign problem. The IKKT model in d = 10 dimensions is also called the IIB matrix 
model. It is postulated to give a constructive definition of type IIB superstring theory. 

Mass deformations such as the Myers term [iQl are essential in order to reproduce 
non-trivial geometrical backgrounds such as the fuzzy sphere in Yang- Mills matrix models. 
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Supersymmetric mass deformations in Yang-Mills matrix models and Yang- Mills quantum 
mechanics models are considered for example in |14lll5j . Yang- Mills quantum mechanics 
models such as the BFSS models [3] in various dimensions are a non-trivial escalation 
over the IKKT models since they involve time. The BMN model [B] which is the unique 
maximally supersymmetric mass deformation of the BFSS model in d = 10 admits the 
fuzzy sphere background as a solution of its equations of motion. The BFSS and BMN 
models are postualted to give a constructive definition of M-theory. 

The central motivation behind these proposals of using Yang-Mills matrix models 
and Yang-Mills quantum mechanics as nonperturbative definitions of M-theory and su- 
perstring theory lies in D-brane physics [150H152] . At low energy the theory on the 
{p + 1)— dimensional world- volume of A^ coincident Dp-branes is the reduction to p + 1 
dimensions of 10 dimensional supersymmetric Yang-Mills [4]. Thus we get a (p+1) dimen- 
sional vector field together with 9 — p normal scalar fields which play the role of position 
coordinates of the coincident N Dp-branes. The case p = corresponds to DO-branes. 
The coordinates become noncommuting matrices. 

As we have already said the class of matrix models of interest to us in this article are 
IKKT matrix models in d = 3 and d = 4 dimensions with mass deformations. The main 
reason behind this interest is that these matrix models suffer generically from an emergent 
geometry transition and as a consequence they are very suited for studying nonperturba- 
tively gauge theory on Moyal-Weyl spaces. Furthermore the supersymmetric versions of 
these matrix models provide a natural nonperturbative regularization of supersymmetry 
which is very interesting in its own right. Also since these matrix models are related to 
large A^ Yang-Mills theory they are of paramount importance to the string/gauge duality 
which would allow us to study nonperturbative aspects of gravity from the gauge side of 
the duality. These motivations can also be found elsewhere, for example in |153ftl56] . 

The motivation for studying large N Yang-Mills matrix models with backgrounds given 
by fuzzy spaces is therefore four-fold: 

• 1) This is the correct way of defining nonperturbatively noncommutative gauge the- 
ories using random matrix models. 

• 2) This provides also a nonperturbative definition of supersymmetry. Indeed super- 
symmetry in this language may prove to be tractable in Monte Carlo simulation. 

• 3) They provide concrete models for emergent geometry. Indeed geometry in tran- 
sition is possible in all these matrix models. 

• 4) These large N Yang-Mills matrix models are also relevant for the string/gauge 
duality. 

Among the tools that can naturally be used in the context of Yang-Mills matrix models are 
Monte Carlo simulation [15711158] . 1/N expansion [132^ll59ltl61j . renormalization group 
equation [162^1163] and random matrix theory |143[|144lll64] . Monte Carlo simulation 
using the Hybrid Monte Carlo algorithm was adapted to this type of matrix models 
in [46 yi65j . However simulations of matrix models remain much harder than simulations 
of field theory because of the non-local character of matrix interactions. Renormalization 
group approach to matrix models can be found for example in |166j and |167ffl.71j . 
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The central result in this review section is as follows. In the case of gauge theory 
on fuzzy complex projective spaces |172j we will describe how the corresponding matrix 
models allow for a new transition to and from a new high temperature phase known as 
Yang-Mills or matrix phase with no background geometrical structure. The low tempera- 
ture phase is a geometrical one with gauge fields fluctuating on a round complex projective 
space. We discuss the case of the fuzzy sphere in great detail then towards the end we 
comment on the other known cases. 

Noncommutative gauge theory on the fuzzy sphere was introduced in [14711148] . As 
we have already mentioned it was derived as the low energy dynamics of open strings 
moving in a background magnetic field with S^ metric in |86H88j . This theory consists of 
the Yang-Mills term YM which can be obtained from the reduction to zero dimensions of 
ordinary U{N) Yang-Mills theory in 3 dimensions and a Chern-Simons term CS due to 
Myers effect |16j . Thus the model contains three N x N hermitian matrices Xi, X2 and 
X3 with an action given by 

5 = YM + CS = -\Tr[Xa, Xi,f + '^eabcTrXaX^^X,. (4.2) 

This model contains beside the usual two dimensional gauge field a scalar fluctuation 
normal to the sphere which can be given by |2UUj 

* = 4^. (4.3) 

The model was studied perturbatively in |59j and in [60 1I173J . In particular in [59] the 
effective action for a non-zero gauge fluctuation was computed at one-loop and shown to 
contain a gauge invariant UV-IR mixing in the large N limit. Indeed the effective action 
in the commutative limit was found to be given by the expression 

+ non local quadratic terms. (4-4) 

The 1 in 1 -|- 2g(^A3 corresponds to the classical action whereas 2g(^A3 is the quantum 
correction. This provides a non-local renormalization of the inverse coupling constant 
1/ g^. The last terms in (j4.4p are new non-local quadratic terms which have no counterpart 
in the classical action. The eigenvalues of the operator A3 are given by 

£^^lGl + l)^2(/2 + l)^ 
X (/2(Z2 + l)-/l(/l + l)))- 




p{p+ 1 



/=1 



In above L + 1 = N. The 1 in 1 — (— iyi+'2+P corresponds to the planar contribution 
whereas (— 1)'i+'2+p corresponds to the non-planar contribution where p is the external 
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momentum. The fact that A3 7^ in the hmit N — > means that we have a UV-IR 
mixing problem. 

The model YM + CS was solved for A^ = 2 and iV = 3 in [174]. It was studied 
nonperturbatively in [58J where the geometry in transition was first observed. 

In |57j a generalized model was proposed and studied in which the normal scalar field 
was suppressed by giving it a quartic potential V with very large mass. This potential on 
its own is an 0(3) random matrix model given by 



V 



N 



2C2 



a^lJiTriXl) 



(4.6) 



The parameter /i is fixed such that // = m . The model S + V was studied in [55] and [56] 
where the instability of the sphere was interpreted along the lines of an emergent geometry 
phenomena. For vanishing potential m^ , ^ — > the transition from/to the fuzzy sphere 
phase was found to have a discontinuity in the internal energy, i.e. a latent heat (figure 
[1]) and a discontinuity in the order parameter (figure [2]) indicating that the transition is 
first order. The order parameter is identified with the radius of the sphere, viz 



_ — TrD X 

r Nc2 



aDa 



(4.7) 



From the other hand the specific heat was found to diverge at the transition point from 
the sphere side while it remains constant from the matrix side (figure [3]) . This indicates 
a second order behaviour with critical fiuctuations only from one side of the transition. 
This to our knowledge is quite novel. The scaling of the coupling constant a in the large 
N limit is found to be given by a = ayN . We get the critical value 



a, = 2.1 ±0.1. 



(4.8) 



The different phases of the model are characterized by 



fuzzy sphere (a > a* 
r = 1 
C, = 1 



matrix phase (d < d*) 

r = 

a = 0.75 



For ?n, 7^ and/or /U 7^ the critical point is replaced by a critical line in the ji — t plane 
where /3^ = d^/(l + m?Y and t = /i(l + m?). In other words for generic values of the 
parameters the matrix phase persists. The effective potential in these cases was computed 
in [59]. We find 



Feff = a^ 



[i*^ 



1 



+ 



m 



^^21 



3" 4 

The extrema of the classical potential occur at 



+ ln( 



1 



1 + m? 



0, 



i±^/rT4t 



y± 



(4.9) 



(4.10) 
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For fj, positive the global minimum is (j)^. The is a local maximum and (/>_ is a local 
minimum. In particular for /i = m^ we obtain the global minimum 0+ = 1. For /i negative 
the global minimum is still (/>+ but becomes a local minimum and cp- a local maximum. 
If /i is sent more negative then the global minimum 0+ = 1 becomes degenerate with 
= at t = — I and the maximum height of the barrier is given by V- = /3'^/324 which 
occurs at !?!)_ = o. The model has a first order transition at t = —2/9 where the classical 
ground states switches from i?i>+ for t > —2/9 to for t < 2/9. 

Let us now consider the effect of quantum fluctuations. The condition V^^ = gives us 
extrema of the model. For large enough a and large enough m and fi it admits two positive 
solutions. The largest solution can be identified with the ground state of the system. It 
will determine the radius of the sphere. The second solution is the local maximum (figure 
HI of Ves and will determine the height of the barrier. As the coupling is decreased these 
two solutions merge and the barrier disappears. This is the critical point of the model. 
For smaller couplings than the critical value a^, the fuzzy sphere solution Da = (pLa no 
longer exists. Therefore the classical transition described above is significantly affected 
by quantum fluctuations. 

The condition when the barrier disappears is V^g = 0. At this point the local minimum 
merges with the local maximum (figure H]) . Solving the two equations V^g = V^g = yield 
the critical value 

2 1 _ (PlicP, + 2fl) 



at 8 



where 



1 + m^ 



i + ,,i + Ml±^ 



(4.12) 



If we take ^ negative we see that g* goes to zero at /u(l + m?) = — 1/4 and the critical 
coupling d=f is sent to infinity and therefore for /i(l + m^) < —j the model has no fuzzy 
sphere phase. However in the region — ^ < fi{l + m'^) < — | the action S + V is completely 
positive. It is therefore not sufficient to consider only the configuration Da = (pLa but 
rather all SU{2) representations must be considered. Furthermore for large a the ground 
state will be dominated by those representations with the smallest Casimir. This means 
that there is no fuzzy sphere solution for /u(l + m^) < — |. 
The limit of interest is the limit /u = m? — )-oo. In this case 



This means that the phase transition is located at a smaller value of the coupling constant 
a as m is increased. In other words the region where the fuzzy sphere is stable is extended 
to lower values of the coupling. 

Nonperturbatively the value ag defined as the value of a at which curves of the average 
value of the action < 5 > for different A'' cross gives a good estimate of the location of 
the transition. For large A^ we observe that the location cimax of the peak in C„ and the 
minimum dmin coincide and agree well with d^. By extrapolating the measured values of 
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Omax and cimin to A^ = oo we obtain the critical value otc- The critical coupling determined 
either as a.c or as dtg gives good agreement with (I4.13p . The phase digaram is given in 
figure m 

In the case when we include a potential term with /i = m^ we also found numerical and 
analytical evidence |175tll76] for the existence of another transition on the fuzzy sphere 
which is of the Gross- Witten type |177j . This is a field theory transition which occurs 
within the fuzzy sphere phase before we reach the matrix phase. We also note that a 
simplified version of our model with V quartic in the matrices, i.e. vr? = and ^ 7^ was 
studied in pTSlfTTO] . 

In |38] an elegant pure matrix model was shown to be equivalent to a gauge theory on 
the fuzzy sphere with a very particular form of the potential which in the large N limit 
leads naturally (at least classically) to a decoupled normal scalar fluctuation. In [175pi76] 
and [39] an alternative model of gauge theory on the fuzzy sphere was proposed in which 
field configurations live in the Grassmannian manifold U{2N)/{U{N + 1) x U{N — 1)). 
In [39] this model was shown to possess the same partition function as commutative gauge 
theory on the ordinary sphere via the application of the powerful localization techniques 

[MIEZ]. 

The matrix phase which is also called the Yang-Mills phase is dominated by commuting 
matrices. It is found that the eigenvalues of the three matrices Xi, X2 and X^ are 
uniformly distributed inside a solid ball in 3 dimensions. This was also observed in higher 
dimensions in |62j . The eigenvalues distribution of a single matrix say X3 can then be 
derived by assuming that the joint eigenvalues distribution of the the three commuting 
matrices Xi, X2 and X3 is uniform. We obtain 

p{x) = ^,{R'-x'). (4.14) 

The parameter R is the radius of the solid ball. We find numerically the value R = 2. 
A one-loop calculation around the background of commuting matrices gives a value in 
agreement with this prediction. 

These eigenvalues may be interpreted as the positions of DO-branes in spacetime fol- 
lowing Witten |3]. In [61] there was an attempt to give this phase a geometrical content 
along the same lines. However our notion of geometry in this article follows Connes [23] 
which requires providing the Dirac or Laplacian operator together with the algebra in 
order to determine the geometry of the space. Therefore for all practical purposes this 
phase has no geometry since we can not identify a sensible Laplacian acting on commuting 
matrices. 

We recognize two different scaling limits. In the fuzzy sphere phase the matrices Xa 
define a round sphere with a radius which scales as A^ in the commutative limit, i.e. they 
define a two-dimensional plane whereas in the matrix phase they define a solid ball in 
3 dimensions. This is the scaling limit of the flat two dimensional plane. In the second 
scaling limit the scaled matrices Da define on the other hand a round sphere with finite 
radius in the fuzzy sphere phase whereas in the matrix phase they give a single point. 

The essential ingredient in producing this transition is the Chern-Simons term in the 
action which is due to the Myers effect. This transition is related to the transition found 
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in hermitian quartic matrix models. For example the 0(3) matrix model given by the 
potential V does not have any transition but when the Chern-Simons term is added to it 
we reproduce the one-cut to the two-cut transition. By adding the Yang-Mills terms, i.e. 
by considering the full model we should then obtain a generalization of the one-cut to the 
two-cut transition. Indeed the matrix to the fuzzy sphere transition is in fact a one-cut 
to N-cut transition (figure [6]). 

The matrix phase should be identified with the one-cut (disordered) phase of the 
quartic hermitian matrix model. By formally comparing (jl.2p and ()4.6p we can make the 
identification a = —jia^ and b = m? /{2c2)- The one-cut phase corresponds to the region 
a > a=i, = —2-\fh or equivalently 6? < a^ where 

a, = -^. (4.15) 

For fj, = TV? this is precisely the critical point (j4.13p . 

By using the methods of cohomological field theories and topological matrix models 
employed in |33p34p68j . it might be possible to bring the model into the form of a hermitian 
matrix model with generalized interaction of the form |203j 

S = NaTrM^ - Tr^ ln(adM + b). (4.16) 

In summary we find for pure gauge models with global S'0(3) symmetry an exotic line 
of discontinuous transitions with a jump in the entropy, characteristic of a 1st order 
transition, yet with divergent critical fiuctuations and a divergent specific heat with critical 
exponent a = 1/2. The low temperature phase (small values of the gauge coupling 
constant) is a geometrical one with gauge fields fluctuating on a round sphere. As the 
temperature increased the sphere evaporates in a transition to a pure matrix phase with 
no background geometrical structure. These models present an appealing picture of a 
geometrical phase emerging as the system cools and suggests a scenario for the emergence 
of geometry in the early universe. 

Lastly we remark on the effect of fermionic determinants on the transition which is the 
subject of the remainder of this article. We propose here mass deformed supersymmet- 
ric matrix Yang- Mills which are reduced from mass deformed supersymmetric Yang-Mills 
quantum mechanics as the prime candidates to study the effect of supersymmetry on 
emergent geometry and vice versa in noncommutative gauge theory. These mass deformed 
matrix models or quantum mechanics provide also the prime examples of supersymmet- 
ric models which can be put on a computer. It is conjectured that supersymmetry will 
remove the transition and stabilizes completely the geometry against the quantum fluctu- 
ations of the noncommutative gauge theory or else supersymmetry may be dynamically 
broken. Indeed our Monte Carlo results reported here (see section 7) confirms this pic- 
ture. In [U] a Monte Carlo simulation of a 4 dimensional model with a Chern-Simons 
term was performed. Although the model was not invariant under (flat) supersymme- 
try transformations the effect of the added Majorana fermions seemed to stabilize the 
geometry. 

Finally we make few remarks on known analogous results in 4 dimensions. In 4 di- 
mensions we have 2 fuzzy projective spaces, fuzzy CP^ and fuzzy S^ x S^. Classical 
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gauge theory on fuzzy CP^ and fuzzy S^ x S^ are considered in [1801I181J and |182yi83j 
respectively. A Monte Carlo study of the model YM + CS on fuzzy CP^ was conducted 
in |184j . It is observed that in this model both a fuzzy CPjy phase and a fuzzy sphere 
phase exist together with the matrix phase. The phase structure is therefore much richer. 
This was confirmed in [181) with the calculation of the one-loop effective potential of 
the model YM + CS + y on fuzzy CP^. The one- loop effective potential of the model 
YM + CS + y on fuzzy S^ x S^ was computed in |185j and a Monte Carlo study of 
the same model but with V = was performed in [186] with similar conclusions. Fuzzy 
S^ X S^ is also considered in p^7Hl89] 
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Figure 1: The observable ^^ for m? = as a function of the coupling constant for different 
matrix sizes N. The solid line corresponds to the theoretical prediction using the local minimum 
of the effective potential. 
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5 Mass Deformation of D = 4 Super Yang-Mills 
Matrix Model 

5.1 Dimenional Reduction in 4D 

We work with the metric t] = { — \- ++). The gamma matrices satisfy {7^, ^u} = '^'Hfj.v 
We consider the representation 



We have /3 = 74 = ^7°, 75 = i7o7i7273- We verify that 7+ = 7°7^7°, 7* = 7^7/^7^, 
7J = C^fiC. The charge conjugation matrix is defined by 

C = 72/3 = -e75 , C+ = C-i = -C , C^ = -C , C^ = -1. (5.2) 



The Majorana condition reads 

^j} = iilj+jO = -ilFC. (5.4) 

The N = 1 supersymmetric Yang-Mihs theory is given by the Lagrangian density 

A = Tr ( - \f,,F^'' - \\rD,\ + \d'^ . (5.5) 

Ft,y = df,A^ - duA^ - ig[A^, A^\ , D^\ = d^X - ig[Af„ A]. (5.6) 

The supersymmetric transformations are explicitly given by 

60 A^ = 67^ A 



SoX = 

60D = iej^j^Df^X. (5.7) 



The reduction of this theory to one dimension is obtained by setting A^ = X^/g where 
diX^ = 0. We also set A = 'ip/g and D = F/g where diX = 0, diF = 0. We get the 
supersymmetric Yang-Mills quantum mechanics given by the Lagrangian density 



^Trl^-{DoX,f + \[X,,X,f - ^V^7°AV' + \H\XrM + \f^ 



Dq = dQ - i[XQ, ]. (5.9) 
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The super symmetric transformations become 

5qF = -ie757oDoV' + eJ5'yi[Xi,^p]. (5.10) 

Note that the variation of the Lagrangian density Cq under the supersymmetry transfor- 
mations is given by (we set F = for simphcity) 

r'-i^Y[X„i^] - Tr^^j'Doi; + ^{DoXif + i[ 
1 



3^ 

This vanishes by Fierz identity. 



e^Trtpa{tpu,tp(3}- (5.11) 



5.2 Deformed Yang-Mills Quantum Mechanics in 4D 

Let /x be a constant mass parameter. A mass deformation of the Lagrangian density 
Co takes the form 

2 

C^ = Co + -^Ci + ^C2 + ... (5.12) 

gZ gZ 

The Lagrangian density Cq has mass dimension 4. The corrections £i and £2 niust have 
mass dimension 3 and 2 respectively. We recall that the Bosonic matrices Xq and Xa have 
mass dimension 1 whereas the Fermionic matrices Vi have mass dimension |. A typical 
term in the Lagrangian densities Ci and C2 will contain Uf Fermion matrices, nj, Boson 
matrices and nt covariant time derivatives. Clearly for Ci we must have |n/ + nb + nt = 3. 
There are only three solutions {nf,nh, nt) = (2, 0, 0), (0, 3, 0), (0, 2, 1). For £2 we must have 
^nj + nb + nt = 2 and we have only one solution {nf,ni),nt) = (0,2,0). Thus the most 
general forms of £1 and £2 are 

£1 = Tr Umx/j + ^SabcXaXbX, + JabXaDoXh) . (5.13) 

C2 = Tr(-^SabXaxX (5.14) 

Clearly for £3 we must have |nj + n^, + nt = 1 which can not be satisfied. Thus the 
correction £3 and all other higher order corrections vanish identically. 

We will follow the method of [15] to determine the exact form of the mass deformation. 
We start with the fermionic mass term 

^Cj, = ^TrijMtl; = -fiTri^CM^Cip. (5.15) 

gZ -V gZ 
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We can verify that only the identity matrix, the gamma five and the cubic terms in the 
gamma matrices can survive in the expansion of M. The cubic terms are either 7*^ [7*, 7-'] 

or 7^7^7^. we have 

C^ = Tri^UaU + 675 + ^i/ij7°[y,7-''] + C7S'7^] V'- (5.16) 

Under the chiral transformation ip — > x = e^'^'^^ip, -ip — > x = 'ipe'^'^^^ the Lagrangian 
density Cq remains the same while C^ transforms as 

£^ = TrxfiaU + fo'75 + ^H^.j^'b',!'] + nW^X, (5.17) 

where a = acos 20+6sin2(/> and b = — asin2(/> + 6cos20. Thus we can use this symmetry 
to set 6 = 0. We get 

£^ = Tri,(iaU + ^^y7°[y,7''] + C7S'7^]V'. (5.18) 

The numerical coefficients a, Hij and c will be constrained further under the requirement 
of supersymmetry invariance. 

Next we consider the bosonic terms. We can choose the coefficients Jab to be antisym- 
metric without any loss of generality since we have f dt TrJabX^-DoX^ = f dt Tr-^iJab — 
Jba)X'^DQX^ — f Tr-^DoJabX^'X where clearly the last term can be included in £2- The 
Bosonic part of the action reads 

l-Trf^DoXaf + l[Xa,Xbf + ^SabcXaXbX, + fiJabXaDoXb - ^SabXaXb] .(5.19) 



(7^ \2 4 3! 2! 

We consider the action of a time dependent 50(3) rotation L defined by Ya = LabXb, 
Xa = LbaYb, LL = L L = 1. The action transforms as 

^Trl^iD^Yaf + ][Ya. Ybf + ^^SabcYaYbY, + fiJabYaDoYb - ^SabYaYb] , (5.20) 
where babe = Jaoboco^aao^bbo^ccoi ^ab = '^aofeo^aao^f'fco ~ u'^«obo^aaoC'0^fe6o ~ 7^^0-^acC'O-t^fec; 

Jab = Jaob^LaaoLbbo + j^LbcdoLac- Thus it is clear that we can choose L such that Jab = 0. 
By similar arguments we can show that the coefficients Sab can be chosen to be totally 
symmetric while the coefficients Sabc can be chosen totally antisymmetric. By rotational 
invariance we must therefore have Sab = v6ab and Sabc = Qie^abc for some numerical 
coefficients v and e. In particular the Myers term is 



— ^myers = — f leCijkTrXiXjXk j . (5.21) 



The mass deformed super symmetric transformations will be taken such that on bosonic 
fields they will coincide with the non deformed supersymmetric transformations so that the 
Fierz identity can still be used. The mass deformed supersymmetric transformations on 
fermionic fields will be different from the non deformed supersymmetric transformations 
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with a time dependent parameter e = e(t) which satisfies Oqc = //He. We will suppose the 
supersymmetric transformations 

^^Xq = 60 Xq 
5^Xi = SoXi 
Sf^ip = 6qiP + ;uAe. (5.22) 

By requiring that the Lagrangian density ()5.12p is invariant under these transformations 
we can determine precisely the form of the mass deformed Lagrangian density and the mass 
deformed supersymmetry transformations. A long calculation yields the mass deformed 
Lagrangian density and mass deformed supersymmetry transformations given respectively 
by (see Appendix C) 

C^ = Co + —Tr'tlj{ia- —-f^j'^j^)ip + ieeijk—TrXiXjXk ^-{e^ - —a^)TrXf. 

(5.23) 



6f,Xo = €70-0 
6^Xi = e-iitp 

e = 6(t) = e'^(^^"+t^V7V)*. (5.25) 

We verify that {6^X^)^ = 5^X^ and hence the Hermitian matrices X^ remains Hermitian 
under supersymmetry. The corresponding supersymmetric algebra is su{2\l) jl5] . 

We prefer to work with the parameters /^i and ^2 defined by ia/x = /Ui/4 and — |^ = 
/i2/4. The Lagrangian density and supersymmetric transformations become 

C^^l = ^o + ^Tri;{fii+H2l^7^7^)il^ -ieijk:^TrXiXjXk- Y^{^^'i + nl)TrXf. 

(5.26) 

Sfj,Xo = e-fotp 
5^Xi = e-fii; 

S,^=(^- ^[7°,y]^o^« + lb\^n[X^,X,] - ^(/ui - ^27SV)y^.)e.(5.27) 



(/xi70-^27°7S^7^)t 



e = e(t) = e6l'^i^'-^27"7 7 7";t_ ^^^S) 
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5.3 Truncation to Zero Dimension 

We consider now the Lagrangian density (action) given by 

Cf, = £o + ^Tr'i[){fii + fj.2j^-f^'j^)ip-ieijk^TrXiXjXk- -^{fil + f4)TrXf. 

(5.29) 



^Tr(^^[X„X,][X'^,Xn+'-' 



Co = -Tr( -[X^,X,][X'^,X'^] + -i,j'^[X^,ij] ] . (5.30) 



In above we have ahowed for the possibihty that mass deformations corresponding to the 
reduction to zero and one dimensions can be different by including different coefficients a, 
b and c in front of the fermionic mass term, the Myers term and the bosonic mass term 
respectively. However we will keep the mass deformed supersymmetric transformations 
unchanged. We have 

5^X^ = ej^ip = -ip^f,e. (5.31) 

S,i^ = (^[7^71[^M>^-]-^(w-M27^V)f^^)e• (5.32) 

We compute the supersymmetric variations 



g^6^Co = -^rrV;(/Ui +;U27S'7')7°7i^o,^.]e- ^rrVS(/^i -/i27'7'7')[y ,7^1[^.,^,]e- 



(5.33) 



<5^(^^TrVi(/xi + /i27W)V') = ^{fJ.l + fiDTrd^iXf 



+ ^rrV;(m + /^2y7V)7V[^o,Xi]e 

+ '-^Tri;{^l^ + ^,2l'7S')h\7'][X^,X,]e. (5.34) 

Clearly the first term of ()5.33p must cancel the second term of ()5.34p . i.e. 

2 
a = - , c = 1. (5.35) 



We get then 



5^ (^g^Co + ^rr^(/ii + fi2l'lW)i^ " l^^^' + f^D^rXf^ 



^- ^' .1^,2^3 Uo,* V 



12 



rr^(/ii + /X27V7')[7^7^■][^^,^i]e• 



(5.36) 
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By using the identity ^7"^7^7^[7*,7-'] = —^ijkl^ we find 

S.I^-ie^.J'-fTrX^X.X^ = -'-^Tri,^^j^j'[f,^^][X„X,]e. (5.37) 

We must then have 

6 = 1. (5.38) 

But also we must have 

/ii = 0. (5.39) 

Thus we get 

5^(^g^£o + ^Tri,j'lh^i^-^TrXf-ie,,k^TrX,X,Xk] = 0. (5.40) 
The model of interest is therefore 



3 



ieijk^TrX.X.Xk). (5.41) 



Since ip and e are Major ana spinors we can rewrite them as 



^4-r •-rrn- 



We compute with Xq = iX/^ the action 



Cf, = ^Tr(l[X^,X,]^ + U[X,,Xj]-Ai,kXk] +e+(i[X4,..]+ai[X,,..] + ^^ 



(5.43) 
The super symmetric transformations are 
6^Xo = i{uj+e - e+uj) 



5^Xi = i{9^aiuj - uj'^Ui 



1 i 

■:j^ijk(^k[Xi,Xj] + -, 



6^9 = i - i(Ti[Xo,Xi] - -€ijkak[Xi,Xj] + -fi2criXi ]u}. (5.44) 



6 Cohomological Approach 
6.1 Super symmetry Transformations 

The J\f = I supersymmetric Yang-Mills theory in four dimensions is given by the 
Lagrangian density 

A = Tr ( - ^F.^F^'^ - Y^^D.X + ^D^) . (6.1) 
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F^,u = d^A, - d,A^ - ig[A^, A,] , D^X = d^X - ig[A^, A]. (6.2) 

The supersyninietric transformations are explicitly given by 

6A^ = €jfj.X 



6D = ie757^^^A. (6.3) 



'' ■ 4 



The Majorana spinors A and e can be written in terms of two dimensional complex spinors 
6 and u as 

We will also write 

1 2i , , 

Af, = -X^ , D = -B. (6.5) 

We will work with Euclidean metric, i.e Xq = iX/^. The reduction of the above theory 
to zero dimension is obtained by setting d^X^ = 0, d^Xa = 0, d^6 = 0, d^O'^ = and 
d^D = 0. We obtain 

S = -^Tr(j[X^,X,f + ea''[X^,e]-2BA. (6.6) 

We can also trivially check that (we set g^ = 1) 
1 



S = --Tr[X^,X,Y-Tr9+li[X4,..]+aa[X,,..]]9 + 2TrB'. (6.7) 

The supersymmetric transformations become 



69 = ia>'''[X'',X'']uj-2Buj 
59 = -iLoa^"'[X^, X"] + 2Bu 

5B = ]-CJa^'[X^',9] + ]-[X^',9]a^'uJ. (6.8) 

Let us note that since we are in Euclidean signature the transformation law of X4 is 
antihermitian rather than hermitian. 

By using a contour shifting argument for the Gaussian integral over B we can rewrite 
the auxiliary field B as 

B = H + ]^[XuX2]. (6.9) 

The matrix H must be taken hermitian. We will also introduce 

0i = m + im , ^2 = Xi + iX2- (6.10) 
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</. = ^(Xs + iX^) , <^ = -^(Xg - iXi). (6.11) 

Let us now compute 

b = ^iJcohom- (O.izj 

5cohom = Tr Ti^^ + if [Xi , X2] + [X, , <P] [X, ,^] + [cf>,^f-^, [</., r]i] - Xi[^, x^] 

- me'nXi,X,]+m[Xi,xA. (6.13) 



In above we have used e^^ = 1 and the result TrFi[B,F2] = —Tr{Fi,F2}B where Fi are 
fermionic matrices and i? is a bosonic matrix. The indices i and j take the values 1 and 
2. As we will see (p and (j) must be taken to be independent. Furthermore (j) must be taken 
antihermitian and (j) hermitian. 

We have four independent real super symmetries generated by the four independent 
grassmannian parameters ^j , pi defined by the equations uji = ^2 + i^i and uj2 = pi + ip2 ■ 
The supercharges Qi = Qir + iQii and Q2 = Q2R + iQ2i are defined such that the 
supersymmetric transformation of any operator = 0/ + iOr is given by 

60 = [Q+u-uj+Q,0] = {50)R + i{50)i, (6.14) 

where 

{dO)R = -2i[Q, Oi] , {60)i = 2i[Q, Or]. (6.15) 

The generator Q is antihermitian and it can be written as 

Q = ^2QlR + ilQlI + P1Q2R + P2Q2I- (6.16) 

For a bosonic field O the imaginary part O7 is zero and thus we obatin {SO)i = 2i[Q, O] = 
SO = —i50. After a long calculation we obtain 

5X1 = 2(^2X1 + 6X2 + Pim + P2V1) 

6X2 = 2{i2X2 - iiXi + P2m - Pim) ■ (6-17) 

54) = -2{pixi + P2X2) 

54> = -2(6^+6^2). (6.18) 

5H = 2U<P,m]-'^Ci[(P,m] + '^pi[^,X2]-2p2[^,xi] + '^[Xi,pim-p2m] 

+ 2[-pir]2-p2Vi,X2]. (6.19) 

For fermion fields we obtain {56i)r = —5r]i, {59i)i = 6r]2, (<J^2)_R = —^X2 and {662)1 = 
6xi- Another long calculation yields 

-^771 = 2Ci[<P,^]-2C2H-2ipi[Xi,^]-2ip2[X2,^] 

SV2 = -2C2[<PA]-2CiH-2ip2[Xi,^] + 2ipi[X2A]. (6.20) 
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-5x2 = -2p2[(l)A]-2piH-2ip2[X^,X2] + 2^i[X^,(l)] + 2C2[X2A] 

6x1 = 2pi[^,^]-2p2H + 2ipi[Xi,X2]-2C2[Xi,(l)] + 2Ci[X2,^]. (6.21) 

We look at the supercharge associated with S,2- We define the exterior derivative d 
on bosons by dB = i[QiR,B] and on fermions by dF = i{Qifi,F}, i.e 5B = 2^2dB and 
5F = 2S,2dF. The corresponding supersymmetric transformations are 

dX, = Xi. (6.22) 

# = , d^ = -r/2. (6.23) 

dH = [cP,7]i]. (6.24) 

dr]i = H , dr]2 = [^A]- (6.25) 

dxi = [(p,X,]. (6.26) 

From these transformation laws we can immediately deduce that for any operator O we 
must have 

d^0 = [4>,0]. (6.27) 

Thus d'^ is a gauge transformation generated by (j) and as a consequence it is nilpotent on 
gauge invariant quantities such as the action. 
Next we compute 

s(-Trx^[Xi,^]] =2C2Tr([X„ct>][X„^]-x^[^,X^] + V2[x^,xA (6.28) 

s(Trm[Xi,X2]] = 2^2Tr(H[Xi,X2] - m^^J[x^,xA (6.29) 

dfrrmH) = 2i2Tr{H^ - 7?i[0,r/i]Y (6.30) 

5(-Tm2[4>,^]\ =2i2Tr{[ct),4>f-1^2[(t>,m]\■ (6.31) 

Hence 

5TrQ = 2e25cohom. (6.32) 

Q = -xAX^, ^] + ??i[^i, ^2] + ViH - m[(^, 0]- (6.33) 

As a consequence 

dTrQ = S'cohom- (6.34) 

Thus we have 

d^TrQ = dS'cohom = 0. (6.35) 

In the above equation we have used the result dS'cohom = 0. Thus d is nilpotent on gauge 
invariant quantities such as Q which are formed from traces. The term cohomology comes 
precisely from the analogy of d with an exterior derivative. 
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6.2 Cohomologically Deformed Supersymmetry 

We consider the deformed action and deformed exterior derivative given by 

'S'def = ^cohom + S = Scohom + eiSl + €282 + ■■■ (6.36) 

ddef = d + eT. (6.37) 

Super symmetric invar iance requires 

ddciSdci = 0. (6.38) 

The fact that d? is equal on gauge invariant quantities, i.e. d^^cohom = d'^Si = leads 
to d^ S'def = 0. We have the identity 

difSdef = 0. (6.39) 

This is equivalent to 

{d, TjScohom + eT^Scohom + e{d, T}S + eT^S = 0. (6.40) 

Thus we must have among other things 

Kr}Seohom = WT}S = 0. (6.41) 

In other words {d, T} generates one of the continuous bosonic symmetries of the action 
Sdei which are gauge transformations and the remaining rotations given by the 5*0(2) 
subgroup of 50(4). Following p3] we choose {d,T} to be the rotation U defined by 

U : Xa — > ieabXb , Xa — > if^abXb- (6.42) 

We have then 

{d,T} = U. (6.43) 

The symmetry T must also satisfy 

T^ = 0. (6.44) 

By following the method of [48] we can determine precisely the form of the correction T 
from the two requirements T^ = and {d, T} = U and also from the assumption that T is 
linear in the fields. A straightforward calculation shows that there are two solutions but 
we will only consider here the one which generates mass terms for all the bosonic fields. 
This is given explicitly by 

TXi = , Txi = ie-ijXj , T(j) = 

TH = i-fr]2 , Tr]2 = , Trji = -iXcf) + i-f^ , T^ = 0. (6.45) 
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The cohomologically deformed super symmetric transformations are therefore given by 

ddeiXi = xt- (6.46) 

ddef'/' = , ddefA = -m- (6-47) 

ddef-ff = [0,'^i] + «e7r/2. (6.48) 

ddeiVi = H + e{-iX(p + i-f(i>) , ddeiV2 = [(!>, (A]- (6-49) 

ddciXi = [(p,Xi]+ieeijXj. (6.50) 

Furthermore we have the result 

dl^^ = cf + eU. (6.51) 

On [/—invariant quantities we have d^^j = (P. For example d^^j = d'^ on quantities 
independent of X^ and Xa- Also d^^f (xa-'^a) = d^{xaXa) and ^^j (e''''xaXb) = d^{e''^XaXb). 
Thus we have on [/—invariant quantities the identity 

dlA-)=d\...) = [<l^,..]- (6.52) 

6.3 Cohomologically Deformed Action 

Next we need to solve the condition (j6.38p . The deformed action is a trace over some 
polynomial P. In the non-deformed case we have S = dQ where Q is a [/- invariant 
expression given by (|6.33p . We assume that the deformed action S^ei = TrP is also 
[/- invariant. By using the theorem of Austing [45] we can conclude that the general 
solution of the condition (I6.38p . or equivalently of the equation ddcfTrP = 0, is 

•Sdef = ddcfTrQdcf + TrR^{4)). (6.53) 

For SU{N) gauge group this result holds as long as the degree of P is less than 2N/3. 
Clearly when the deformation is sent to zero ddef — > d, Qdei — > Q and R — ;■ 0. Thus 
we take 

Qdei = Q-iR , R = KiRi + K2-R2- (6.54) 

We choose Ri and R2 to be the [/- invariant quantities given by 

^1 = ^f^abXaXb , R2 = -Vi4>- (6-55) 

We choose R3{<j)) to be the [/- invariant quantity given by 

Rsi^) = -P^<P^- (6.56) 
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In order to remove the deformation we must take e — > so that djcf — > d and p — > 
so that S'def — > dTrQdci and kj — > so that Qdcf — > Q- 
We compute 



S'def = dTrQ- idTrR + eTTrQ - ieTTrR + TrR^ {(j)) 

— '-'cohoin I "J • 



(6.57) 



The first term Scohom = dTrQ is the original action. By using the result that d{FiF2) 
dFi.F2 - FidF2 and T{FiF2) = TF1.F2 - F1TF2 we obtain 



idTrR = —iTr 



Ki(l)[Xi,X2] - K.1X1X2 - K2H(j) - K2r]ir]2 



(6.58) 



eTTrQ = eTr 



i\ct>{H + [X^,X2]) + n^H - i77?ir/2 + i(7 + 2)^[Xi,X2] 



. (6.59) 



ieTTrR = eTr I -niXf + K2X(p4> — K.2l<P' 



Also 



TrR3{(p) = Tr{-p 



2j,2\ 



(6.60) 



(6.61) 



We will choose the parameters so that the total action enjoys 5*0(3) covariance with 
a Myers (Chern-Simons) term and mass terms for all the bosonic and fermionic matrices. 
The relevant fermionic terms are 



iKiTr{xiX2 - mm) + i{i^i + H2- e-i)Trr]ir]2 = iKiTr{xiX2 - mm)- 
For SO{'i) covariance we have chosen 

/^i + '^2 - £7 = 0. 
Next we remark 

Tr(H'^ + H[Xi,X2]\ + i{K2 + e-f)TrH^ - ieXTrH(P = 

Tr(^H + ^[Xi,X2]+'-iK2 + e^)^-'-eX^^ - ^Tr[X,,X2f + 

^Tr ({k2 + e-/)^- eX(Pj -^Tr({K2 + £7)^ " exA [X^ , X2] . 

The relevant bosonic mass terms are therefore given by 

Tr i -eniXf + eK2A(/)(^ - eK2^(iP - p^cj)^ I + -Tr ( {k2 + e7)<A - '^X4> j 

-eKiTrXf + (-e^A^ - p^)Tr4>'^ + -(ks - e-ffTr^'^ + -eA(K2 - ej)Tr(P^. 
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(6.62) 



(6.63) 



(6.64) 



(6.65) 



In order to cancel the cross product we choose 



-/X'-p' = \{n,-e^f. (6.66) 



Then 



Trl -€KiXf + eK2X4>4' - e«^27</'^ - /'^'/'^ ) + j^ri (k2 + (^7)4> - eA(/) 

-eKiTrXf + -{k2- e^){K2 - ej - e\)TrXl - -{k2 - e-/){K2 - e-y + eXjTrXJ = 
z o o 

^eKiTrX^ -1{k2- ej){K2 - ej + e\)TrXl (6.67) 
z o 

Again for 5*0(3) covariance we have chosen 

-e^i = -(k2 -e7)(K2 -e7-eA). (6.68) 

The Chern-Simons term can be obtained from the terms 
Tr(- iiKi + eA)<^[Xi,X2] + ie(7 + 2)^[Xi,X2]] - '-Tr ({^2 + e^)^ - eX^ [Xi, X2] = 

-'-{eX + 2Ki)Tr(l)[Xi,X2] + ^{-K2 + ej + 4e)Tr^[X^,X2] = 
- J(eA + e7 + 2ki - K2 + 4e)TrX3[Xi, X2] + -(eA - 67 + 2ki + K2 - 4e)TrX4Xi,X2].{6m) 
The first term is the Chern-Simons action. We want to impose the condition 

eA + 2ki + k2 = £7 + 4e. (6.70) 

The solution of equations IK63h . ^Mh . ^Mh and ([670]) is 

j = ^l±^ , A = 4-^, p2 = 2e(2e-^i). (6.71) 

e e 

The deformation action AS'cohom is therefore 

1 1 

AS'cohom = «'ti7'r(xiX2 -??ir/2) + -e^iTrX^ - -(k2 -e7)(K;2 -e7 + eA)rrX| 

z O 

- ^{eX + e-f + 2Ki-K2 + 4e)TrX3[Xi,X2] 

= '^Tr0+0 + ^e^iTrX^ + 1^1(26 - Ki)TrXl - '-{Ae + Ki)eafecTrX,XfeXe. 

(6.72) 

We introduce now 

(X 

- (4e + Ki) = a , Ki = -- + 4Coa. (6.73) 



Thus we get 



« /. .„ ^ x„ „J^^ o 



A5eohom = -7(1 - 12Co)Tr0+0 + ^(1 + 6Co)(l - l2Co)TrX', + — Co(l - 12Co)TrXi 
D 3d z 

+ ^aeafceTrXaXftXe. (6.74) 

6 
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The total action is 



Sdei = (Scohon. - Tr{H^ + H[Xi,X2])] +Tr(H + ^[Xi,X2] + ^{k2 + 67)^ - '-exA 

- ^Tr[Xi,X2f + AScohom. (6.75) 

The first hne is effectively equivalent to Scohom- Thus the total action is 

'-'def = '-'cohom + '^'-'cohom- (0.7dJ 

Next we perform the scaling 



X^ -^ {2N)-^X, , e -^ J^^^e, (6.77) 

and 

a — > 2{2N)^a. (6.78) 

We get then the one-parameter family of actions given by (we set B = 0) 



-TT\X^,X,f + N —ea^^TrXaX^X, + __ ^^ > -^ ^/1 -^ ^^.^^2 



5dcf = —rTT\X^, X,f + iV^-e,b,rrX,X,X, + ——(1 + 6Co)(l - 12Co)rrX^ 



+ 4iVa2Co(l - 12Co)TrX| - -Lrr0+ ^[^4, ..] + aa[Xa, ..] + ^(1 - 12Co)V- 

(6.79) 
For stability the parameter (^0 must be in the range 

< Co < 1/12. (6.80) 

This action for Co = is precisely the mass deformed action derived in section 2. The 
value Co = 1/12 will also be of interest to us in this article. This one-parameter family of 
actions preserves only half of the A/" = 1 super symmetry in the sense that we can construct 
only two mass deformed supercharges [38| . 



7 Simulation Results for i^ = 4 Yang-Mills Matrix 
Models 

7.1 Models, Supersymmetry and Fuzzy Sphere 

We are interested in the cohomologically deformed Yang-Mills matrix models 
- i[X^, X^f + ^e,,,X,XfeXe + N^TrXl + Np^TrXJ 
-Tre+(i[Xi, ..] + aa[Xa, ..] + c] 0. (7.1) 



NTr 

1 

Wa 
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The range of the parameters is 



/? = ^(a + 6eo)(a - 12^0) , P^ = 4Co(a - 12^0) , C = ^(« - 12^0). 



(7.2) 



o<eo< 



a 
12' 



(7.3) 



This action preserves two supercharges compared to the four supercharges of the original 
non deformed Yang-Mills matrix model [IB]. We will be mainly interested in the "mini- 
mally" deformed Yang-Mills matrix model corresponding to the value ^o = a/12 for which 
we have 



5 



NTr 



T [^ti 1 ^y] H :r' ^abcXaXbXc 



1 

Wa 



Tre+{i[Xi,..]+aa[Xa,.. 



(7.4) 



The "maximally" deformed Yang-Mills matrix model corresponding to the value ^o = 
coincides precisely with the mass-deformed model in D = 4 and as such it has a full 
Af = 1 mass deformed supersymmetry besides the half M = 1 cohomologically deformed 
supersymmetry. From this perspective this case is far more important than the previous 
one. However there is the issue of the convergence of the partition function which we will 
discuss shortly. In any case the "maximally" deformed Yang-Mills matrix model is given 
by the action 



S 



NTr 
1 



1 2 2ia 

T [Xfj, , X^\ H — eabcXaXhXc 



+ n'^-^TtXI 



Na 



Tr9+{i[X^,..]+aa[Xa,..] + -a 



(7.5) 



The above two actions can also be rewritten as (with Xa = aDa and a = ayN) 



'S'suSY 



NTr 
1 



T [Xa , Xi)] H — eabcXaX^Xc 



+ NJ3a^TrXl 



Na 



Tr9+{i[X^,..]+aa[Xa,..]+a^ 



(7.6) 



Equivalently 



SsVSY 



a^ 

N' 
1 

N^ 



Tr 



^[Da,Dbf + jeabcDaDbD, 



1 ^"^'r n2 
+ j;fTrDa 



Tre-^[i[Di,..]+aa[Da,..]+i 



(7.7) 



/3 = , ^ = cohomologically deformed 



/3 



e 



mass deformed. 



(7.8) 
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We remark that the bosonic part of the mass-deformed Yang-Mihs matrix action can be 
rewritten as a complete square, viz 



Sb = NTr 



T [^t^ ) ^u] H :j- ^abcXaXhXc 



+ N^-^TrXl 
9 



= NTr(^[X^,X,] + ^e^,^X^ . (7.9) 

Clearly e^^x = if any of the indices /U,z^,A takes the value 4. Generically the bosonic 
action of interest is given by 



Sb = NTr 



jK..^J + f 



7 [X^ , X,^] -\ — EabcXaXhXc 



+ Ni3a^TrXl. (7.10) 



Here we allow (3 to take on any value. The variation of the bosonic action for generic 
values of /3 reads 

5Sb = NTrJ^SXi + NTrJ^dXh 

Ja = [Xa, [Xa, X^]] , J, = 2^a^X, + i[Fab, Xa] + [X4, [X4, X,]] , 

Fab = i[Xa,Xb]+aeabcXc. (7.11) 

Thus extrema of the model are given by 1) reducible representations J^ of SU{2), i.e 
Xa = Ja and X4 = and 2) commuting matrices, i.e X^ belong to the Cartan sub- 
algebra of SU{N). The identity matrix corresponds to an uncoupled mode and thus we 
have SU{N) instead of U{N). Global minima are given by irreducible representations of 
SU{2) of dimensions N and 0. Indeed we find that the configurations Xa = (pLa, X4 = 
solve the equations of motion with (j) satisfying the cubic equation 0(<^^ — a<p + /3a^) = 0. 
We get the solutions 



1±a/1-4/? 
(/>o = 0, 0± = a "^ . (7.12) 

We can immediately see that we must have /?<l/4 which does indeed hold for the values 
of interest /? = and /? = 2/9. However the action at (f)± is given by 

S*J = ^o^(a-iT^;^^). (7.13) 

We can verify that 5b [(^_] is always positive while Sb[4>+] is negative for the values of 
(3 such that /?<2/9. Furthermore we note that S[(/>o] = 0. In other words for /3<2/9 
the global minimum of the model is the irreducible representation of SU{2) of maximum 
dimension A'^ whereas for /3 > 2/9 the global minimum of the model is the irreducible 
representation of SU{2) of minimum dimension 0. 

At /3 = 2/9 we get 0+ = 2a/3 and S'[0+] = 0. Thus the configuration Xa = ^La 
becomes degenerate with the configuration Xa = 0. However there is an entire SU{N) 
manifold of configurations Xa = ^ULaU~^ which are equivalent to the fuzzy sphere 
configuration. In other words the fuzzy sphere configuration is still favored although now 
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due to entropy. Thus there is a first order transition at /3 = 2/9 when the classical ground 
state switches from Xa = -^ La to Xa = as we increase /3 through the critical value 
(3 = 2/9. The two values of interest /3 = and /? = 2/9 both lie in the regime where the 
fuzzy sphere is the stable classical ground state. 

This discussion holds also for the full bosonic model in which we include a mass term 
for the matrix X4. Quantum correction, i.e. the inclusion of fermions, are expected to 
alter significantly this picture. 

Towards the commutative limit we rewrite the action, the second line of (j7.7p . into the 
form (with a^ = 1/g'^ and Fab = i[Da,Db] + €abc4>Dc) 

^ = Wn^'^'' " W^' V^^'^P^^D, + '^^a] + ^ (</'(<A - 1) + /?) TrDl 
- ^^^[^«'^4]' - ^Tre+{i[D^,.]+<Ja[Da^.]+i)0. (7.14) 

The 3rd terms actually cancel for all values of /3. Thus 

S = -l^TrFl-^-^Tr[eabcFabD. + ct>Df^-^Tr[Da.D,f 

- ^Tre^{i[D^,..] + aa[Da,.]+i)0- (7.15) 

The commutative limit A^ — >• 00 is then obvious. We write Da = (piLa + Aa) and we 
obtain 



S = ^lj!.-^^^eab.lFabA.-^ljCaD,f-l^^ 

(7.16) 

7.2 Path Integral, Convergence and Observables 

In the quantum theory we will integrate over N x N bosonic matrices X^ and N x N 
fermionic matrices 6*+ and 9a- The trace parts of X^, 9^ and 9a will be removed since they 
correspond to free degrees of freedom. The partition function of the model is therefore 
given by 

Z = f VX^ V9 V9^ 6{TrX^) 6{Tr9^) 6{Tr9a) e'-^susY 

= I VX^6 {TrX^) detV e'^''. (7.17) 

detP = I d9d9+5{Tr9a)5{Tr9t,)e^'^''^^^^. (7.18) 

The integration over the fermions yielded the determinant of the 2{N'^ — 1) x 2(A^^ — 1) 
dimensional matrix T) = i[X4, ..] + aa[Xa, ■■] + ^a. This determinant is positive definite 
since every eigenvalue A of P is doubly degenerate [46] . The reason lies in the fact that 
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the Dirac operator V = 1X4^ — iX^ + (TaXa — CTa^a + C'^ i^ symmetric under the exchange 
of left and right operators, viz under Xa -H- —X^. A much cleaner proof goes as follows. 
Let ^ be an eigenstate of T> with eigenvalue A, in other words 

i[X4,^] + aa[Xa,^]+ia7p = X^. (7.19) 

Taking the hermitian conjugate of this equation we get 

^[^4, i^+f] - a'^iX,, i^+f] + ia{^+f = X{^+f. (7.20) 

In above (^+) is a column vector with components given by '^^2- Multiplying the above 
equation by (T2 and defining the spinor ^ = cr2(^^) we arrive at the equation 

i[Xi,^] + aa[Xa,^]+Ca^ = X^. (7.21) 

We have also used the identity aa = —'720" J 0-2. We conclude that ^ is also an eigenstate 
of T> with the same eigenvalue A. The spinors ^ and ^ are charge conjugate to each 
other. In above we have assumed that A is real which follows from the fact that the Dirac 
operator D = i[X4, ..] + aa[Xa, ..] + ^a is hermitian. This establishes that the determinant 
defD is positive definite for any configuration X^ and hence the model can be accessed 
directly by Monte Carlo simulation. 

Let us also add that the Dirac operator D admits an approximate chirality operator 
and hence there is an approximate chiral symmetry in this model beside exact rotational 
invariance, exact gauge invariance and exact charge conjugation. The existence of chiral 
symmetry (though approximate) means that there should exist more structure in the low 
energy fermionic spectrum. 

The partition function Z is also invariant under the translation X^ — > X^ + eX^ 
where e is a small parameter. Under this coordinate transformation the measure dX^ 
changes to (1 + 4(iV^ — l)e)(iX^. The bosonic action Sb = S4 + S3 + S2 changes to 
Sb + e(4S'4 + 353 + 2S'2) under X^ — > X^ + eX^. The determinant, on the other hand, 
changes under X^ — > X^ + eX^ as 

detP -^ {l + ef^^'-^'^det(M' -eia{l+^)] 

-^ (l + e)2(^'-i)detYl--^ea(l + 7))det'A^'. (7.22) 



The matrices A4 and 7 are given in appendix B. We obtain then 
detP ^ ( 1 + e 2(iV2 - 1) - iaTr'^^^ - iaTr'^^^j 



deW. (7.23) 



From the invariance of the partition function under the coordinate transformation X^ — > 
X^ + eX^ we derive therefore the Schwinger-Dyson identity 

,<YM> ,<YMo> ,<CS> ^<RAD> -COND 6 

IDE = 4 K h4 ^ h3 ^ \-2 ^ \-fa ^ \ ^ 

= 6. (7.24) 
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This is an exact result. 

The operators YM, YMq and CS are the actions given by 

YM = -^Tr[Xa,X,f , YMo = -^Tr[Xo,Xaf , CS = 'J^eabcTr[Xa,X,]X,. 

(7.25) 

The action RAD is related to the radius of the sphere. It is given by 

RAD = Npa^TrXl (7.26) 

We will define the radius r of the sphere through the relation 

The total bosonic action is given by S" = YM + YMq + CS + RAD. The condensation is 
defined by 

COND = -^^ =< ^Tre+e > . (7.28) 

We find 

COND = CONDo + CONDi. (7.29) 



CONDo =< Tr^d — >> 

CONDi =< Trld ^ 7 > • (7.30) 

4X4,. .J ^Oa{Xa,.\ +4" 

These are real quantities. The condensation CONDq is obviously real since every eigen- 
value of the Dirac operator is doubly degenerate. The expression of the total condensation 
COND suggests that it is real. We have verified numerically that CONDi is real and that 
the condensation is dominated by CONDq. 

For completeness we will also measure the logarithm of the determinant of the Dirac 
operator as an independent observable in Monte Carlo simulation. 

Convergence of Yang- Mills path inetgrals such as the one given by ()7.17p was studied 
extensively in [15] and in |49fl54] . This question is of paramount importance for analytical 
analysis as well as for Monte Carlo simulation. The source of the divergence, if any lies 
in the so-called flat directions, i.e. the set of commuting matrices. 

The path integral (j7.17p corresponds to a gauge theory with gauge group SU{N) in 
dimension D = 4. We will also consider SU{N) gauge theory in dimension D = 3. We 
start the discussion with the model a = 0, /3 = and ^ = 0. It was found in [l9] that the 
bosonic path integral in D = 3 is convergent for A^ > 4 while the bosonic path integral in 
D = 4 is convergent for A^ > 3. Since we are interested in large values of A^ we can safely 
consider the bosonic path integrals in Z) = 3, 4 to be convergent for all practical purposes. 
On the other hand, it is found in ^53j, that the supersymmetric path integral in D = 3 
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is not convergent while the supersynimetric path integral in D = 4 is convergent for all 
iV> 2. 

Tuning the parameters a, /3 does not change this picture. For example it was shown 
in [S3] that adding a Myers term, i.e. considering a non-zero value of a, does not change 
the convergence properties of the D— dimensional Yang-Mills matrix path integral. The 
point is that the Chern-Simons (Myers) term is always small compared to the quartic 
Yang-Mills term. The same argument should then lead to the conclusion that adding a 
bosonic mass term, i.e. if we consider a non-zero /3, will not change the above picture. 

Tuning the fermion mass term, i.e. considering a non-zero value of the scalar curvature 
^, will lead to complications. In this case the Pfaffian, or equivalently the determinant, 
will be expanded as a polynomial in the scalar curvature ^. The analysis of [53] should 
then be repeated for every term in this expansion. We claim that the super symmetric 
path integral in D = 4 is not convergent for generic values of ^. 

We have extensively checked in Monte Carlo simulation the conjecture that Yang-Mills 
matrix models in dimension D = A does not make sense for generic values of ^. The major 
observation is that for i_^ the fermion determinant for generic values of a = \/Na never 
reaches thermalization [j. However, we have also observed that for sufficiently small values 
of ^ the theory actually makes sense and thus there is some critical value of ^, which we 
will not determine in this article, above which the path integral is ill defined. The value 
of interest ^ = 2/3 corresponding to the mass deformed matrix model lies in this range 
where the model is actually ill defined. 

Therefore in order to access the mass deformed Yang- Mills matrix model by the Monte 
Carlo method we must regularize the theory in such a way as to make sure that the path 
integral is absolutely convergent. Unfortunately most regularizations will not maintain 
neither the full M = 1 mass deformed supersymmetry nor the half A/" = 1 cohomologically 
deformed supersymmetry of this model. We adopt here the regularization in which we 
will simply set ^^ = 0. In other words we make the replacement 



1 2 2za 



+ N^a^TrXl 
Tre+[i[Xi,.]+aa[Xa,.]]- (7-31) 



'S'SUSY — > 'S'suSY = XTr 

1 

In summary, the value j3 = 2/9 corresponds to the mass deformed Yang-Mills matrix model 
with softly broken supersymmetry whereas the value /? = is precisely the minimally 
deformed model which enjoys half of the N = 1 cohomologically deformed supersymmetry. 
The Metropolis algorithm and other algorithms used to study these models numerically 
and the detail of the simulations can be found in appendices A and B. 



^This happens typically for small values of a far from the fuzzy sphere region but not very close to a = 0. 
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7.3 Bosonic Theory: Emergent Geometry and Phase Dia- 
gram 

Emergent Geometry We measure the different observables as a function of the cou- 
pUng constant a for the two relevant values of /3, i.e. /3 = 0, 2/9. We have verified that 
the bosonic Schwinger-Dyson equation holds in Monte Carlo simulations. See figure ()12p . 
Recall that this identity reads for the bosonic models as follows 

TDF - . <YM> <YMo> ^ < CS > < RAD > 4 _ 

The radius which we have defined by the equation < 1/r >=< TrX\ > /cP'C2 is shown 
on figure ([7]). For large values of a the result is consistent with the classical prediction 



1 TrX"^ o l + Vl-4/3 
< - >=< -^ >= <^^ , 0+ = h^ -• (7.33) 

This means in particular that the system is in the ground state configurations 

X^ = , Xa = a(t>La. (7.34) 

In other words we have a fuzzy spherical geometry given by the commutation relations 

[X4,Xa] = , [Xa,Xb] = ieabca(pX,. (7.35) 

We have checked these commutation relations and found them to hold quite well for 
sufficiently large values of a. The coordinates on the sphere are defined by 

na = ^ ,^nl = (^\ (7.36) 



C2a 

We observe that as we decrease a, the radius 1/r jumps abruptly to then starts to 
increase again until it becomes infinite at d = 0. This is the most dramatic effect of the 
so-called sphere-to-matrix transition in which the sphere suddenly expands and evaporates 
at the transition points then it starts shrinking to zero rapidly as we lower the coupling 
further. 

This is the interpretation advocated in |55H57j for a similar phenomena observed in 
the case of three dimensional bosonic models. As far as we know this phenomena was 
observed in Monte Carlo simulation first in [58] and it was found in analytical work on 
perturbative three dimensional bosonic models in [59] and then in [60| . 

The transitions for the bosonic mass deformed model with /3 = 2/9 and the bosonic 
cohomological model with (3 = are observed to occur at the following estimated values 



a^ = 4.9 ±0.1 , /3 = 2/9. (7.37) 

a* = 2.55 ±0.1 , /3 = 0. (7.38) 
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The fuzzy sphere phase corresponds to the region a > a^, whereas the matrix phase 
corresponds to the region d < q*. In other words the sphere becomes more stable as we 
make (3 smaher (see below). 

The order of the sphere-to-matrix transition is very difficult to determine. Since the 
ground state configurations are X4 = and Xa = a(l)La, the theoretical analysis based 
on the effective potential of the three dimensional model done in |55y56j should also hold 
here largely unchanged (see below). As a consequence we will only summarize here the 
main points omitting much technical details. 

The specific heat C^, = (<5^> — <S' >^)/(A^^ — 1) shown on figure ([8]) diverges 
from the side of the fuzzy sphere with a critical exponent equal 1/2. It is equal to 3 in this 
phase where 1/2 is due to the 2 dimensional U{1) gauge field on the sphere, 1/2 is due to 
the normal scalar field on the sphere and 1/2 is due to the scalar field X4. This critical 
behavior is typical of a second order transition. In the matrix phase the specific heat is a 
constant right up to the transition point and it is equal 1 where each matrix contributes 
1/4. There is no divergence from this side and the critical exponent is 0. In other words 
the behavior above and below the critical coupling are different, which is quite unusual, 
but still from the specific heat we qualify this transition as second order. 

The expectation values of the Yang-Mills action and the Myers (or Chern-Simons) 
action are shown on (|10p . The expectation values of the total bosonic action for the two 
cases f3 = 2/9 and /3 = are shown on figure (jlip . From these observables we observe a 
discontinuity at the transition point. Thus the sphere-to-matrix transition is associated 
with a latent heat equal toA<S'>=<5 >matrix — < S >sphere which is typical of a 
first order phase transition. It is straightforward to estimate the value of this latent heat. 
The latent heat is released by going from the matrix phase to the fuzzy sphere phase for 
/3 = whereas for /3 = 2/9 the latent heat is released by going in the other direction from 
the fuzzy sphere phase to the matrix phase. 

As we will see from the discussion of the eigenvalues distributions the matrices X^ in 
the matrix phase are commuting matrices centered around 0. 

Phase Diagram The last point we would like to address within the context of the 
bosonic theory is the construction of the phase diagram in the plane a — 13. We have 
already measured two points of this phase diagram which correspond to the two values 
/3 = 2/9 and /3 = 0. In order to map the phase boundary between the fuzzy sphere and 
the matrix phase we choose other values of /3 and measure for each one of them the critical 
value of a from the discontinuity in the radius 1/r. The result is shown on figure ([9]). 

The effective potential of the 4 dimensional bosonic Yang-Mills matrix model in the 
Feynman-'t Hooft background field gauge in the ground state configurations X4 = and 
Xa = a(j)La can be calculated using the method of [59]. We find 



VcE ~4 

= a 

2C2 



4 3 2 



+ 21og</.2. (7.39) 



The difference with the three dimensional bosonic Yang-Mills matrix model lies in the 
factor of 2 multiplying the logarithm. The critical line can then be obtained following the 
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method of [55]. We get 



^. = 3 ^ L _ 5^) . ^J = ^i^_. (7,40) 

This is the fit used on figure ([9]) with very reasonable agreement with the Monte Carlo 
data. 

7.4 Dynamical Fermions: Impact of Supersymmetry 

In this section we will discuss the effect of the fermion determinant. First we note that 
simulations with fermions are much more harder than pure bosonic simulations. The main 
source of complication is the evaluation of the determinant which is highly non trivial. 
Thermalization is very difficult and as a consequence taking the limit of large A^ is not 
so easy even with the use of the Hybrid Monte Carlo algorithm. In the bosonic case we 
could go as large as A^ = 100 with very decent number of statistics although in this article 
we have only reported data with A^ up to A^ = 16. In the fermionic case we will report 
data with A^ up to A^ = 10. 

The first thing we have checked is the Schwinger-Dyson identity (j7.24p . The Monte 
Carlo data are shown on figure ()12p . The fermionic data agrees well with the prediction 
6 whereas the bosonic data agrees well with the prediction 4. Note that 6 = 4 + 2 where 
4 is the number of bosonic matrices and 2 is the number of fermionic matrices. 

The most important order parameter with direct significance to the underlying geom- 
etry is the radius 1/r. See figures ()13p . (I14p and (J15l) . We observe that the transition 
sphere-to-matrix observed in the bosonic theory disappeared completely. Again it seems 
here that there is no major difference between the two models with f3 = 2/9 and /3 = 0. 
It is clear from the structure of the action that the theory with a = can not support 
the fuzzy sphere geometry and thus a phase with commuting matrices must still exist. 
However, the transition to the phase of commuting matrices starting from the fuzzy sphere 
phase seems to be a crossover transition not the second/first order behavior observer in 
the bosonic model. This seems to be confirmed by the behavior of the specific heat, the 
Yang-Mills and Myers actions and the total action shown on figures ()16p . p7til8p and ()19p 
respectively. The jump and critical behavior in the specific heat and the discontinuity in 
the various actions disappeared. 

We have to note here that the observable < TrX^ > /N diverges in the supersymmetric 
theory with a = /3 = [51j. For the mass deformed theory we have /3 = 2/9 and thus the 
observable < TrX"^ > /N always exists. We observe on the second graph of figure (I13p 
that < TrX^ > /N increases as we decrease a towards which is consistent with the fact 
that it will diverge in the limit a — > 0. Qualitatively the same phenomena is observed 
for /3 = on the second graph of figure (|14p with more erratic behavior as we decrease 
a towards 0. However in this case we can not infer that < TrX"^ > /N exists for all a 
since /3 = although it looks that it does from the data. From this perspective the mass 
deformed model is better than the cohomologically deformed model. 

We have not succeeded in determining precisely the value at which the crossover tran- 
sition occurs but it seems that it depends on N in such a way that it is pushed to smaller 
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values of a as we increase A^. From figure (J15p we can see that the crossover transition for 
/3 = 2/9 occurs at d = 3.13,2.63,2.38 and 2.13 for A^ = 4, 6,8 and TV = 10 respectively. 
A simple fit yields the result 

"* = ^ ' /3 = 2/9. (7.41) 

The conjecture that the crossover transition occurs at arbitrarily small values of a in the 
large A^ limit is one of the main results of this article. In this way the fuzzy sphere is truly 
stable in the super symmetric theory and does not decay. In any case we are certain that 
the fuzzy sphere in the supersymmetric theory is more stable compared to the bosonic 
theory and the crossover transition to the matrix phase is much slower. This conclusion 
is similar to that of [47j . 

7.5 Eigenvalues Distributions 

Bosonic Theory A powerful set of order parameters is given by the eigenvalues dis- 
tributions of the two matrices X3 and X4. See (PU|) . The eigenvalues distribution of the 
matrix X4 is qualitatively the same for all values of the coupling constants d. However, 
the eigenvalues distribution of the matrix X3 suffers a major change as we go across the 
transition point. In the fuzzy sphere region the eigenvalues distribution of X3 is given 
by an A^— cut distribution corresponding to the A'^ eigenvalues — (A^ — l)/2, ..., (A^ — l)/2 
whereas in the matrix phase the eigenvalues distribution of ^3 is identical to the eigen- 
values distribution of A'4. 

The eigenvalues distribution p/^{xi^ of the matrix X4 is always centered around 0. In 
the fuzzy sphere phase p/^{x/^ depends on the coupling constant a. In the matrix phase 
below the critical value the eigenvalues distribution pdjyx^ does not depend on a and 
coincides with the eigenvalues distribution of the non deformed model with d = 0. In this 
region the eigenvalues distributions of the matrices X3 and X4 are identical. 

Motivated by the work |6H I62j it was conjectured |201j that the joint eigenvalues 
distribution of d matrices X\, X2^...Xd with dynamics given by a reduced Yang-Mills 
action should be uniform inside a solid ball of some radius R. We have already checked 
that this conjecture works in three dimensions [63]. We will check now that this conjecture 
holds also true in four dimensions. Let /o(a;i, X2, X3, X4) be the joint eigenvalues distribution 
of the 4 matrices Xi, X2, ^"3 and X^. We assume that p(a;i,X2,X3,X4) is uniform inside 
a four dimensional ball of radius R. The normalization condition gives p(x\^X2^xj,^x/^ = 
I/V4 = 2/it'^R^. We want to compute the eigenvalues distribution of a single matrix, say 
X4, which is induced by integrating out the other three matrices. We compute 
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We obtain therefore the eigenvalues distribution 



(7.42) 



P4ix4.) 



3ttR^ 



{R^ 



Xa 2. 



(7.43) 



This is precisely the fit shown on figure ([20]) with a measured value of i? for /3 = and 
a = given by 



Ro = 1.826 ±0.004. 



(7.44) 



The above eigenvalues distribution works better for the theory with /3 = 2/9 as shown on 
figure ()22p with a similar measured value of R given for a = 0.25 by 



R. 



■2/9 



1.815 ±0.008. 



(7.45) 



We have found that these two measured values are almost the same throughout the matrix 
phase. 

We emphasize that P4,{x4^) is the eigenvalues distribution of the matrix X4 not only in 
the matrix phase but also in the fuzzy sphere phase with a value of R which depends on 
the coupling constant a. We also emphasize that Piix^) is the eigenvalues distribution 
of the matrix X3 in the matrix phase for the bosonic theory for both values /3 = and 
/3 = 2/9. 

Another non trivial check for this important conjecture is the theoretical prediction of 
the radius 



< ^TrX! > 



3 < ^TrX^ > 

3 / dX4P4^{x4)x 
J-R 

R^ 



(7.46) 



This means that in the matrix phase the order parameter < j^TrX^ >, which is related 
to the radius, is constant. Indeed, this is what we see on figure ([7]) with a good agreement 
between the Monte Carlo measurement and the theoretical prediction. The observed value 
of < jjTrX^ > is slightly below the theoretical prediction for /3 = 2/9 throughout. This 
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is not the case for (3 = where the Monte Carlo measurement starts shghtly below R'^/2 
and then rises above it as we approach the transition to the fuzzy sphere. 

This can potentially be a serious difference between the tow cases /3 = 2/9 and (3 = 0. 
The transition to the sphere in the case of (3 = 2/9 is in the form of an abrupt jump but 
in the case of /3 = there is a slow rise in the matrix phase as we increase a before the 
actual jump. 

The main conclusion of these successful measurements is the fact that the matrices 
Xfj_ in the matrix phase are commuting and thus they are diagonalizable with a joint 
eigenvalues distribution which is uniform inside a ball of dimension R. 

Supersymmetric Theory In the supersymmetric case we found it much easier to 
compute eigenvalues distributions with the value (3 = 2/9 and thus we will only consider 
here the mass deformed model. A sample of the eigenvalues distributions of the mass 
deformed model is shown on figures (j2ip and (|22p . 

Again it is observed that the eigenvalues distribution of X3 in the fuzzy sphere phase is 
given by an A''— cut distribution corresponding to the A^ eigenvalues —{N — l)/2, ..., (A^ — 
l)/2 whereas in the matrix phase the eigenvalues distributions of A3 is given by p^^x^) 
with a much larger value of R given for a = 0.25 by 

ii = 2.851 ± 0.009. (7.47) 

The eigenvalues distribution of X4 is always centered around given by ^4(2:4) with a value 
of R which depends on the coupling constant a. This eigenvalues distribution coincides 
with the eigenvalues distribution of A3 deep inside the matrix phase below a = 0.25. 

The eigenvalues distribution p^ is therefore universal in the sense that it describes the 
behavior of the eigenvalues of the matrix A4 for all values of a and all values of (3 and 
the behavior of the eigenvalues of the matrix A3 in the matrix phase for all values of (3. 
We note here the difference between this eigenvalues distribution p4 and the eigenvalues 
distribution of the supersymmetric model with a = (3 = [51]. In the latter case the 
distribution extends from —00 to +00 and goes as 1/x^ for large eigenvalues. It is not 
clear to us at this stage how the two distributions relate to each other. 

The point at which the eigenvalues distributions of A3 and A4 coincide may be taken 
as the measure for the crossover transition point and thus for A^ = 10 this occurs at a 
point between d = 1 and d = 0.25. 

7.6 Remarks: D = 3 Yang-Mills Matrix Models and Scalar 
Fluctuations 

The D = 3 Yang-Mills matrix models we can immediately consider here can be ob- 
tained from the above D = A models by simply setting the fourth matrix A4 to zero. This 
is different from the IKKT supersymmetric Yang-Mills matrix model in D = 3 by the 
fact that it involves a determinant instead of a Pffafian and as a consequence the path 
integrals of these theories are convergent. 
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The physics of these D = 3 models is identical to the physics of the D = 4 models in 
the sense that there is a first/second order phase transition from a background geometry 
(the fuzzy sphere) to commuting matrices which in the presence of dynamical fermions 
is turned into a slow crossover transition. The most important difference is the natural 
expectation that the eigenvalues distributions of the matrices X^ in the matrix phase must 
be distributed according to the formula 

P3ix3) = ^iR'-xl). (7.48) 

By analogy with the L) = 4 formula (j7.43p this distribution can be derived from the 
assumption that the joint eigenvalues distribution of the 3 matrices Xi, X2, X3 is uni- 
form inside a three dimensional ball of radius R. In the next section we will also derive 
this distribution for the D = 4 Yang-Mills matrix model with a particular choice of the 
parameters with (3^ ^ (see (j7.ip ). 

Monte Carlo measurment of the radius gives the value 



i? ~ 2. (7.49) 

As shown on the first graph of the figure (j23p a sample of the data for the A'^ = 10 three- 
dimensional bosonic Yang-Mills matrix model with /? = 2/9 and a = 0.5 is shown. Clearly 
it can be fit nicely to (|7.48p with R = 2. In performing the fitting in three dimensions 
we have to cut the tails in order to get a sensible answer. As it turns out the same three 
dimensional data can also be fit to the four dimensional prediction (I7.43p . We also show 
the data for the A = 10 four-dimensional bosonic Yang-Mills matrix model with /3 = 2/9 
and d = 0.5 for comparison. 

In the Monte Carlo data of the D = 3 Yang-Mills matrix models reported in this 
article it was not possible to resolve the ambiguity between the two fits ()7.43p and (I7.48p . 
However high precision runs performed in |202| seems to indicate that indeed the three 
dimensional prediction (j7.48p is the correct behavior for the eigenvalues distributions of 
the D = 3 Yang-Mills matrix models. 

The second remark we would like to discuss in this section concerns the dependence on 
A^ and a of the eigenvalues distributions /04 and ps given in ()7.43p and (|7.48p respectively. 
As shown on the second graph of the figure (|23l) the distributions /94 and p^ are independent 
of a. Similarly we can show that these distributions are independent of A. 

The third remark concerns the eigenvalues distributions of the normal scalar field in 
the fuzzy sphere phase which is define by (/> = (A^ — (p'^C2)/{2(py/c2) |59j . This is shown on 
figure ()24p for both the D = 3 and D = 4 models. The behavior in both cases is the same 
although we have to note that the effective values C2^eS in Z) = 3 and D = 4 are slightly 
different. This is quite natural as the three dimensional model is more stable in the sense 
that it has a lower critical value d=K . The central observation here is that we can nicely fit 
these eigenvalues distributions to the eigenvalues distribution p^ given in (|7.43p as shown 
on figure (IM]) . 
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8 Elements of the Solution of the D = 3 Theory 
8.1 Localization: Reduction to D = 3 Chern-Simons Theory 

We start by summarizing our main results of section 3. The action of the theory 
of interest consists of two terms Scohom and S, viz 5def = -^cohom + >§ u- The first piece 
5cohom is the standard M = 1 Yang- Mills matrix theory in D = 4 which has 4 supercharges 
whereas the second piece S is the cohomological deformation which preserves two of the 
supercharges. The parameters of the deformation were fixed so that the action (after the 
integration of the BRST field H) corresponds to a generalization of the mass deformed 
action considered in section 2 with mass terms for the bosonic and fermionic matrices 
and a Myers term. Effectively the action 5'dcf will be the sum of two actions 5cohom and 
A5cohom given by 

5cohom = -^Tr[X^,X,f-^Tre+(i[Xi,..] + aa[Xa,..]]e + TrB\ (8.1) 

A^cohom = ^Tr9+9 + -eKiTrXl + -Ki{2e - KiJTrXJ - ^(4e + Ki)eabcTrXaXhX^. 
2 2 4 



The first supercharge of this cohomologically deformed theory corresponds to the super- 
symmetry transformations 

ddeiXi = Xi- (8.3) 

ddef0 = , ddef0 = -r]2- (8-4) 

ddeiH = [(/',r/i] +i(K2 +Ki)r/2. (8.5) 

ddcim = H + (-i(4e - ki)(/) + i{K2 + ki)0) , d^eim = [(^, </>]• (8.6) 

ddeiXi = [4>,Xi]+ieeijXj. (8.7) 

The second supercharge can be obtained by an appropriate permutation of the spinors r]i 
and Xi PI - 

We observe that we have three parameters ki, K2 and e where K2 does not appear in 
the action. Thus the path integral does not depend on K2- We note that in [48] a different 
cohomological deformation was considered for which it was pssible to show explicitly 
that the path integral is independent of K2. In summary we have then two deformation 
parameters ki and e or equivalently a and Co defined in equation (|6.73|) . For our purposes 



^In fact the action that should be considered is 2Sdef/g^- 
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a will play the role of the gauge coupling constant while Co is the actual parameter of the 
cohomological/mass deformation. 

The theory with K2 = will yield the same action (I8.ip + (l8.2p with the supersymmetry 
transformations 

ddcfXi = x^■ (8.8) 

ddef0 = , (idef<^ = -m- (8-9) 

ddef-H" = [0,?/l] +^'^1^2- (8.10) 

ddef m = -f^ + (-"^(46 - Ki)(/< + iKi(^) , ddcf f?2 = [0, 0] • (8.11) 

ddefXi = [(/>,-'^i] +«eeij-'^i. (8.12) 

We note that the action for the value k\ = 2e corresponds precisely to the mass de- 
formed action considered in section 2 and as a consequence the above supersymmetry 
transformations (I8.8p - (I8.12|) will correspond to one of the mass deformed supercharges. 
The path integral we wish to evaluate is 

Z = j dX^ dB dO^ d0 exp(-5cohom - AScohom). (8.13) 

This can be obtained from the path integral (with S'dcf [^2 = 0] = S'cohom + S[k2 = 0]) 

Z = dXi dX2 dcj) d(j) dH dxi dx2 drji dr]2 exp(— 5dof[K2 = 0]). (8.14) 

In this equation the actions 5cohom and S must be expressed in terms of the BRST fields 
H, (p, (j), Xi, r]i and Xi, viz 

Scohom = Tr{H^ + H[Xi,X2] + [X,,cl)][Xi,^] + [<P,^f -r]i[c|),r]i]-x^[^,X^] 

- rJl^^[xi,X,]+r,2[x^,X,]). (8.15) 

S[k2 = 0] = iKiTr{xiX2-mV2)-^ieTr4>[Xi,X2]+i{Ki + 2e)Tr^[Xi,X2] 

+ iKiTr4>H - i(4e - Ki)Tr(f)H - 2e(2e - Ki)Tr(f)'^ + -cKiTrXf. (8.16) 

To perform the above path integral (j8.14p we will now add the terms proportional to K2 
and since the resulting action remains independent of K2 we will take the limit K2 — > 00 
in which the path integral localizes. We have 

Z = dX^ dB dO^ d9 exp(-5'cohom - AS'cohom) 

= / dXi dX2 d(j) dcf) dH dxi dx2 drji dr]2 exp(— 5def[K2 = 0] — K2AS). (8.17) 
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AS" = iTr7]i7]2 + iTr(l)H + {Ae - Ki)Tr(l)(l) - KiTr(t? . (8.18) 

In the limit K2 — > 00 we are able to integrate the BRST quartet ryi, r/2, H and <p. The 
ferniionic part of the action which depends on K2 is 

iK2Trr]ir]2. (8.19) 

Thus by using the saddle point method we obtain 

dm dm f{vuV2) e— ^^^''1''^ = 4'~h^'-'Nf{0,0). (8.20) 

The actions S'cohom and S[k2 = 0] reduce to 

^cohom = Tr(H^+H[Xl,X2] + [X„4>][X^,^] + [^,^]^-x^[^,Xi]\ (8.21) 

S[k2 = 0] = iKiTrxiX2 - MeTr(j)[Xi,X2] + i(Ki + 2e)Tr^[Xi,X2] 

+ iKiTr^H - i(4e - Ki)Tr<j)H - 2e(2e - Ki)Tr(t>'^ + -eKiTrXf. (8.22) 

Next we do the path integral over H then X^ in that order using again the saddle point 
method. The relevant terms are the bosonic contributions which are proportional to K2. 
They are given by 

(4e — Ki)K2Tr(f)(/) — KiK2Tr(j) + iK2Tr(j)H. (8.23) 

For e < we can verify that the integral over X3 is exponentially damped and therefore 
we can shift X3 appropriately. The resulting integral over H turns out also to be damped 
exponentially. Explicitly we have 

— — — ? O 

{Ae — K,i)K2Tr(j)(j) — KiK2Tr(j)'^ + iK2Tr(j)H = —eK2Tr{^X^ + —{H — KiX^)^ 

- ^Tr{H + (4e - Ki)Xif + ^(4e - Ki)TrXl 
Ibe 2 

(8.24) 

In the limit K2 — ^ 00 the hermitian matrix H will be localized around — (4e — KijXi^ = 
i(4e — Ki)(0 + 0). The above equation reduces to 

{4e — Ki)K2Tr(j)4> — KiK2Tr(j)'^ + iK2Tr(j)H = —eK2Tr{X3 — 1X4)'^ -i (4e — «;i)TrX|. 

(8.25) 

Equivalently 

{Ae — Ki)K2Tr(f)(f) — K,iK2Tr(f) + iK2Tr(f)H = —■ien2Tr(f) {Ae — ni)Tr{(j) + cf)) . 

(8.26) 

We can now shift the integral over X-^ as X^ — > cf) = — ^(^3 — 1X4). Clearly we can 
now assume that (j) is hermitian. For consistency we will also shift the integral over X4 as 
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X4 — > (j) = 1X4^ — (j). In the limit K2 — > 00 the hermitian matrix <j) is locahzed aromid 0. 
The matrix cj) is then seen to be antihermitian identified with 1X4. 
We get thus by using the saddle point method the result 

/ d(h dH fid) H) Qi^^-'^l)'^2Tr^(l)-KiK2Tr4>'^+iK2Tr4>H _ 

4-^^(2.)^^-il/(0,^(4e - Ki)<A)e^(^^--)^^<^l (8.27) 

The K2 dependence cancels completely if we choose 4e = ki. The actions S'cohom and 
S[k2 = 0] reduce to 

Scohom = Tr(-{4e-KifcP^ + i{Ae-Ki)(t>[Xi,X2]\ (8.28) 

S[k2 = 0] = iKiTrxiX2-^ieTr^[Xi,X2] + {4.e-KifTr(t? + 2e{Ki-2e)Tr(t>'^ 

+ ]^eKiTrXl (8.29) 

The integration over the last fermionic degrees of freedom xi and X2 is now trivial since 
they are free degrees of freedom decoupled from everything else. We end up with the 
model 

Z = I dcf) dXi dX2 exp (iKiTr(l)[Xi,X2] - 2e{Ki - 2e)Tr<j? - -eKiTrXf 



+ ^(4e-Ki)rr0M. (8.30) 

This is essentially the path integral of two-dimensional gauge theory on the fuzzy sphere 
studied in [70j. As was shown in [71j it can also be derived from the reduction to zero 
dimension of Chern-Simons theory on S^ . 

We can therefore conclude that the mass deformation considered here reduces the 
4-dimensional supersymmetric theory to a bosonic D = 3 matrix theory which when 
expanded about the ground state yields a noncommutative gauge theory on the fuzzy 
sphere. 

8.2 Exact Integration: Emergence of the Fuzzy Sphere 

In the previous section we have found that (j) = ~{^3 ~ iX4)/2 is localized around 
which means that the most important configurations satisfy X^ — 1X4 = 0. Since 
(j) = 1X4 — (j) we can think of (p as 1X4, i.e. (p is antihermitian, and as a consequence the 
cubic term is not real and the mass term of (p is of the wrong sign. 

Formally we can also think of as a hermitian matrix which can be identified with 
X3. The cubic term in this case is real which is precisely a Chern-Simon term and the 
mass term of (j) is of the correct sign. This can be seen by rewriting in the limit K2 — > 00 
the path integration over X3 and X4 as (see equation (]8.25p ) 

/ dXs j dX4 exp (eK2Tr{X^ - 1X4? " Y(4e - i^i)TrXl\ (...) ~ 

I dX^ l'dX46{Xs - 1X4) exp ( - y(4e - Ki)TrXn (•••)• (8-31) 
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Instead of doing the integral over X3 first then over X4 by performing the shift X3 — > 
(j) = — (X3 — iXii)/2, X/j^ — > (j) = iX/j^ — we will do the integral over X4 first then over 
X3 by performing the shift X4 — > (j) = — (X3 — iX4)/2, X^ — > cj) = X3 + (j). 

In the following we will therefore assume that (f) is a hermitian matrix identified with 
X^. Furthermore recall that e < and as a consequence we will take ki < and ki — 2e < 
in order to have mass terms for (p and Xi. Since e = — a(l + 6Co)/6 and ki — 2e = 6Coa 
we conclude that we must have —1/6 < Co < 0. We have already found that we must 
have < Co ^ 1/12 in order to have a stable classical theory in the direction X4 and that 
the value Co = corresponds precisely to the mass deformed theory of interest which is 
constructed in section 2. 

In the current case we have to take the value 4e = ki, or equivalently Co = ~l/24, 
which is the value at which the dependence on K2 drops completely. For this value the 
field X4 will have a negative mass but the original theory remains stable in the X4 as one 
might easily check in the fuzzy sphere region. 

We introduce 

4 

t = 2e3(Ki - 2e) = -^Co(l + 6Co)'. (8.32) 

io 

Furthermore since (p is localized around we can identify 4> with X3. In order to make 
5*0(3) covariance manifest we perform the rescaling 



The path integral becomes (with 4e = ki) 



Xi -^ -2e\P-^Xi , </> -^ -e(t>. (8.33) 



Z= f d(p dXi dX2 expt(- 2iTr<p[Xi , X2] - Tr<p^ - TrxA . (8.34) 

We perform now the scaling ()6.78p . We get then 

Z = f d(f) dXi dX2 exp Nt( - 2iTr(f)[Xi,X2] - Trcf^ - Trxf\ ■ (8.35) 



t = -^Co(l + 6Co)'. (8.36) 



For 4e = ki we have explicitly 



a^ 



t = — . (8.37) 

32 ^ ^ 



The classical equations of motion are 

[X,,X2] = i^ , [X2, 4>] = iXi , [(/>, Xi] = iX2. (8.38) 

The solutions are given by SU{2) irreducible representations, namely 

Xi = Li, = L3. (8.39) 
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We perform now the integrals over Xi and X2. By integrating over X2 we get 



dXi exp Nt ( - Tr[(t), Xif - Trcj)'^ - TrXf 



Next we diagonalize the hermitian matrix (p. We have 

N 



(8.40) 



1=1 i<j 



iAl) 



We get then 

N 



Z = ffld<p,l[{4>i-(Pjf JdXiexp(-^{Xi)i,{Xi)*j 

4=1 i<j ^ i,j 



Nt{(l)i - <Pjf + Nt 



Nt^, 



Integrating over Xi we get 



142) 



N 
i=l i<j 



•-^i H^j , 



fdQi[ - Nt{(j)i - (j)j)' 



+ Nt] exp(-NtY,4>l\ 



iA3) 



Let us compute the determinant for A^ = 3. Clearly t — > when a — > 0. Thus we have 
in the hmit a — >0 (with Uj = -Nt{(j)i - (pjf + Nt) 



.44) 



det - Nt{4>i - ^j)' + Nt] = 2ti2ti3t23 + in - ^11(^2 + ifs + *2; 



23^ 



But ill is proportional to a. In other words in vanishes if a — > 0. In this case the 
detreminant is given by 



det ( - Nt{(l)i - (pj)^ + Ntj oc JJ ( - iVt 



"J Yj 



cpj)"^ + Nt 



?.45) 



The generalization of this result to higher N is straightforward. This result is actually 
exact and it does not require taking the limit a — ?■ [67j . The path integral becomes 



Z = jj]^d(t>^\{{<Pi-<Pj)^\{(-Nt{(P,-4>,] 

j=l i<j i<j 

This can also put in the form (with t = 1/g^) 

N 



+ Nt\ exp ( - TVt ^ 



'!, — V^J I ^ V ^^ 2 

1=1 J^J -" I 



n#n 



N 



.46) 



^.47) 



68 



We will use the Cauchy formula 



N 



E(-i)°n 



o-eSA 



f{<t>r 



f=l \^i~ ^''{i) +1 



N 






j=i 



+ 1 



(8.48) 



It is quite illuminating to illustrate this identity for smaller values of N . For example we 
consider A^ = 3. We have 






1 



+ 
+ 



H - <Al + 1 02 - 02 + 1 03 - •/'S + 1 

1 1 1 

01 - 03 + 1 02 - 01 + 1 03 - 02 + 1 
1 1 1 

01 - 02 + 1 02 - 03 + 1 03 - 01 + 1 
1 1 1 

01 - 01 + 1 02 - 03 + 1 03 - 02 + 1 
1 1 1 

01 - 03 + 1 02 - 02 + 1 03 - 01 + 1 
1 1 1 



We introduce the notation Ti 



«j 



•"1 



+ 1 



. Then we have 



+ 1 03 - ,/.3 + 1 



/(0l)/(02)/(03). 

(8.49) 



E(-i)"n 



1=1 



/(0i 



K{i) + 1 



/(0l)/(02)/(03) 



(Ti2 + l)(r2i + l)(Ti3 + l)(r3i + 1)(T23 + 1)(T32 + 1) 

- (rf2 - l)(T^l'3 - 1)(T|3 - 1) + (-^13 + l)(Tl2 + l)(r23 + I, 
+ (-Ti2 + l)(-r23 + l)(Ti3 + 1) - {T^2 - l)m'3 - 1) 

- (rf2-i)(r|3-i)-(T23-i)(T|3-i) 

/(0l)/(02)/(03) 

(ri2 + l)(T2i + l)(ri3 + l)(r3i + 1)(T23 + 1)(T32 + 1) 

+ (Ti2 - Ti3 + r23)2 



(T-i2 T~i2 n^^ 
-^ 12^13-^23 



/(0l)/(02)/(03)n 



T, 



V 



i¥=j 



T,j + 1 



^.50) 



By using the Cauchy formula the partition function becomes 

2 = E(-ir/n 



e ^^ 



' 0i - 4'a{i) + 1 



(8.51) 



crg5jv «=1 

From our earlier considerations it is natural to understand the integrals as contour inte- 
grals. The poles are on the real line so we regularize this partition function as follows 



E(-i)"/n 

re5w i=l 



(Te5jv 



- 0a(i) + 1 + i/3 



, /3>0. 



^.52) 
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Thus we close the contours in the upper half-plane. The reason behind this way of 
regularization is twofold. Firstly among the two poles which may appear for each variable 
only one will be counted. Secondly the result obtained here will be consistent with the 
result for Yang- Mills quantum mechanics obtained in |68j. 

Since the eigenvalues (j)i are defined up to a permutation we must also incorporate a 
combinatorial factor equal 1/A^!. Furthermore the tracelessness condition Tr(j) = must 
be included in the form N6{(j)i -|- ... -|- (J)n)- In summary we get the partition function 

Among the A^! integrals there are (N — 1)\ which are non zero. Let us check this in the case 
of iV = 3. The contour integrals over 0i, 02 and 03 of the last three terms respectively 
in (j8.49p vanish because there are no poles in these variables. The first term in (j8.49p 
vanishes because the contour integrals over 0i, (p2 and 03 vanish separately. The second 
and third term in (j8.49p lead to identical contributions given by 





z^ 


= (- 


' 2! 


J5{<t>i 


+ 0S 


+ </'3)il(# 


In general we get 












Zp 


-i 


-1)^- 


{N- 


-1) 
-1) 


/ 


5(01 


N 
+ ... + 0^)J__^ 



e 9 



-i + l + il3 



e 9 



^.54) 



')N- 



(t)i - 0i_i + l + i(5 , 

(8.55) 

Next we will perform the integrals over the variables 02,.--)0Af-i using the residue theorem. 
The two remaining integrals over 0i and 0^^ will be constrained such that 0i -|- 0iv = 0. 
To see this explicitly we try one more time a small value of A^ say A = 4. We have with 
the notation i'j = 1 + i/S the partition function 

^/3 = (-1)^-1 j> 5{(t>i + 02 + 03 + 04) n H'^» ^^_^^_^ + i^ )- (8-56) 

The pole in the 02— plane is at 02 = 03-1-^7 and the pole in the 03— plane is at 03 = cpi+ij. 
By using the residue theorem we get 

2vri)2 /#id04c^(0i +304 + 3.7) e^i^(*?+(^^+^^^)'+(*^+^^)'+<^^) 



4-1 



(-1) 



„J/3(9ii+(04+2i7)+{04+n)+04) . 



1 1 



+ ^7 01 — 04 — 3i7 



>i 



^.57) 



In general there will be a pole in the 02— plane at 02 = 03 + ^7, a pole in the 03— plane 
at 03 = 04 + ij, a pole in the 04— plane at 04 = 05 + i7...and a pole in the 0Af_i— plane 
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at 4>N-i = (f^N + *7- By using the residue theorem we get 

\Ar-l 



-1 



1 



i#7V<5(0i + {N- l)<i)N + -(iV - l)(iV - 2)z7) 



X e 

X 



-■f7(<?^?+(<7^iv+(^~2)i7)2 + ...+(</,jv+i7)^+0^)gj/3(0i+(</,jv+(Af-2)J7)+...+{0jv+i7)+0iv) 
1 1 



4>i- (pN + il (pi - (pN - {N - l)i7_ 

The integration over (pi will be dominated by the pole at (pi 
function becomes 



(8.58) 
N + {N- l)i-f. The delta 



+ {N - 1)(Pn + -{N - 1){N - 2)i-/) = d{N{(PN + -{N-l)i-f) 



In other words we must have (we take the limit /3 



0) 



^.59) 



1 



-2m) 



N-2 



1 



>ia'(PN(- 



.-^+6^) \^-^{^lH^NHN-2)f+...H't>N+i?+^l) 



(Pi-^N + l<Pi-<pN-{N-l) 



(-27ri)^-i 






N 
7 

e 9 



{<P%+{<1>N+(N-2)f+...+{'pN+l?+<1>%) }^ 

N 



N 



9 



E 



m 



JV-l 
2 



iV-1 



= (-l)A^(2vr^)^-ii,exp(-^.(^ + i; 

The smallest eigenvalue is (pN = —{N — l)/2 and the largest eigenvalue is (pi = {N — l)/2. 
We observe that (pi = (pN + A^ — 1. We have in total A'" = 2s + 1 eigenvalues between (p]\[ 
and (pi with a step equal 1, viz m = {N-l)/2, {N-3)/2, ..., -{N-3)/2, -{N-l)/2. This 
vacuum configuration corresponds precisely to the SU{2) irreducible representation s = 
{N — l)/2. The vacuum energy defined as the logarithm of the partition function is clearly 
proportional to the quadratic Casimir in the irreducible representation s = {N — l)/2. In 
summary we have found that the partition function is dominated by the integration in the 



vicinity of the poles 
of SU{2) of size N. 



+ 1 = which corresponds to the irreducible representation 



(8.60) 



8.3 Saddle-Point Method: The Matrix Phase 

Hermitian Case: We will keep assuming that the matrix (p is hermitian. The effective 
potential derived from the path integral (I8.46|) is given by 
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The saddle point associated with this potential is essentially the inverted oscillator problem 
which is the analytic continuation of the super symmetric model considered in [32]. More 
importantly the observation that the saddle point associated with this potential actually 
also corresponds to the Baxter's three-colorings problem [69] . 

As before we will define t = l/g'^. The saddle-point equation reads explicitly 
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N 
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k-(t>i N 
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We introduce the density of eigenvalues and the resolvent given by 



p{x) 
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The saddle-point equation becomes 



dy 



dy 



le 



p{y) 



p{y) 



1 + z -y I- z + y 



+ W{-l-z) + W{l 
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5" J ^ z-y 

In terms of the resolvent this reads 

^z = W(z + ie) + W(z 

r 

We observe that the term X], 1/(1 — (/"fc- 

1/2 whereas the term ^ ■ 1/(1 + (^jt + 

— 1/2. We are interested in the case where the eigenvalues live on a single interval [a, h\ C 

[—1/2, 1/2]. Thus VK(2:) is analytic everywhere in the complex plane except along the cut 

[a, 6]. The above equation can be rewritten as 



-4>j) is singular if one of the eigenvalues approaches 
i)j) is singular if one of the eigenvalues approaches 



W{z + ie) = ^z- W{z - ie) - W{-1 - z) - W{1 - z) , z G [a, 6]. 
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Thus as the complex number z crosses the cut [a, b] from the first sheet into the second 
sheet we see that the resolvent VF(2:) becomes a linear combination of W{z)^ W{\ — z) 
and W^(— 1 — z). In other words we have three cuts [o, 6], [1 — &, 1 — a] and [—1 — 6, —1 — a] 
in the second sheet. By crossing the cuts [a, 6], [1 — 6, 1 — a] and [—1 — b,—l — a] into the 
third sheet we generate more cuts. Hence the domain of definition of Ty(z) is a Riemann 
surface of infinite genus with an infinite number of cuts in each sheet. 

By comparison the one-matrix model (we set 2z/g'^ to V {z)/g^ and M^(— 1 — z) + 
W{1 — z) to 0) yields a Riemann surface with one cut in each of the two possible sheets 
whereas the 0(n)-matrix model (we set Izj g^ to Y {z)/ g^ and VF(— 1 — z) + W{1 — z) to 
nW{—z)) yields a Riemann surface with two cuts in each of the two possible sheets. 

We introduce now the notation z^ = {—l)^{z — k). The saddle-point equation (j8.65p 
is equivalent to the loop equation [67j 



S{z) = 0. 
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f{z) = W^{z) + -(W(l-z) + W(-l-z))W(z)-\(2zW(z)-R(z)). (8.69) 
2 g^ 

The function S{z) has no cut throughout the complex plane and it satisfies S{z+2) = S{z), 
S{1 — z) = S{z). In order to guarantee convergence of the sum we have subtracted the 
polynomial part R{zk) of 2zj.W{zk). 

Next we can easily verify that the resolvent can be rewritten in terms of the moments 
tn of the matrix (j) as 



oo 



Y^ ^n j_ I _,____n _/__\ l 



^(^) = Z^-£^^^-= dxxy{x) = -Trr. (8.70) 

n=0 '' 

After a long calculation we can rewrite the saddle-point equation ()8.67p as 

g— 1 oo oo 

p=0 p=0 m=0 ^ 

It is not difficult to show now that the free theory g^ = Q corresponds to a density of 
eigenvalues given by the Wigner's semicircle law 

p[x) = -^Va2-x2 , a^ = 2g2. (8.72) 

A perturbative solution around this free solution can be constructed following the method 
of |67) . However the exact solution as we will see shortly is quite different from the 
prediction of perturbation theory. This result is rather expected since in the limit of 
interest a — > we have t — > and as a consequence g"^ — > oo. This conclusion is 
actually also confirmed By Monte Carlo simulations. 

Eigenvalues Distribution: As we have already said the analytic continuation of our 
model is essentially the model studied in [34] . However the exact solution we would like to 
construct here is entirely based on the eigenvalues distribution of the matrix (j) from which 
we can calculate all observables of interest to us. This explicit solution is rather different 
from the sophisticated implicit solution found in [34] . The eigenvalues distribution we will 
derive shortly is also different from the one presented in [61] and it agrees well with the 
numerical Monte Carlo simulations. 
Let us first recall our notation 

^"32- 

Let us next consider the saddle point equation once again. This equation can be rewritten 
in the form 

' - ^dy- ^-?M ^, [dyp{y) = l. (8.74) 



t = — = —. (8.73) 



[z-y){l-{z-yf) 
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z /■^+-^ p{z 



^2 
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We have 

dx 1/11 



x{l-x^) 2\ x2 2x4 



+ ... . (8.76) 



dx-^ = ^ + -^3+... (8.77) 

1 — X^ X dX"^ 

h^W^) = 4ln(.^-l). (8,78) 

We will assume 2 things: 

• A quadratic distribution. 

• An (infinite) support [—L,L]. 

We get then to leading contribution in 1/L the expansion 

^ = p(z)(-^ + ...)-p'(z)(| + ...)+/(z)(-| + ...) + ... (8.79) 

We assume a quadratic distribution which is also symmetric given by 

piz) = a + bz'^. (8.80) 

Normalization gives 

1 6L2 , , 

"=21-- '^■*'l) 

The saddle-point equation to leading order in 1/L becomes 

z 6bz 



In other words 

The distribution becomes 
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We imposed the condition 
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By performing the rescaling z — )• zq = {Lq/L)z we can bring the above distribution to 
the form 

po(^o) = ^{Ll-zl). (8.86) 

This can be derived from the same saddle-point equation (j8.74p and therefore the same 
path integral (|8.46p with a coupling constant Iq = 1/ Qq where qq is defined by the equation 

Lo = {\)-*^o- (8.87) 

The original normalization of the matrices corresponds to Lq/L = acj). Thus we have 

Lq = a(j)L 
= 2\/30 

= V3(l + ^1-4/3). (8.88) 

This is independent of N and a which is precisely what we observe in Monte Carlo 
simulations. For the value of (3 at hand which for Qq = —1/24 is /3 = 1/4 we get the 
prediction Lq = v3- If we simply apply this formula for (5 = 2/9 we obtain the prediction 
Lq = 4/-v/3 = 2.31 which should be compared with the measured value Lq = 2.85. 

Since the above distribution will work only for large L it must only be valid for large 
g^ which is equivalent to small a. This is the regime of the matrix phase. We must have 
then 

L»l^g'^»-. (8.89) 

In other words we must have the following lower estimate of the critical value 

a « a* = 2^3. (8.90) 

Antihermitian Case: The above distribution is the same distribution which was 
found for large values of g^ for the antihermitian model in [61]. However the crucial 
difference is the functional dependence of L on (7 which we will now briefly discuss. 

To get the model of [M| and [61] we must assume that (j) is antihermitian and take 
Ki — 2e > and ki < 0, i.e. Co > 0. We get the path integral 

Z = I d(f> dXi dX2 exp Nt( - 2iTr(l)[Xi,X2] + Trcj? - Trxf\ ■ (8.91) 



In this case 



Integration over Xi and X2 yields 

N 



32' 



.92) 



4=1 i<j i<j 



Z = ff[d(t>^ll{(l)i-4>jfll(-Nt{4>i-^jf + Nt\ expTiVt^^ 

(8.93) 
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We get the effective potential (shifting the eigenvalues as (pi — > icpi) 

-V,s{4>i) = -NtY,<Pl + \Y.H^i-<Pj?-lY.^^U<Pi-4>jf + t\{8M) 



The saddle-point equation is 



1 _ 1 ^ 1 



2 - I '^y77—X7T-77—A2-v (8-96) 



Therefore in this case the saddle-point equation is given by 

±= [dv ^ 

g^ J ^(z-y)(l + (z-y)2) 

We go now through the same steps. We have 

z /•^+^ , p(z - x) 

— ' dx- 
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"^ 1 '. . f'^^ , 1 ". . f'^^ , X 
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,_L X{l + X^) ^''J,-L l + X^ ''Jz-L 2(1 + X2) 

lp{z)ln{l + ^)Cl - P'i^) arctanx|^+i + -, 

(8.97) 



-p{z) ln(l + -)|^+i - p {z) arctan xj^+i + -p (z) ln(l + x'Wzt^E + ••. 



Again we assume a symmetric quadratic distribution. As before only the above three 
terms will contribute. The first term is still of order 1/L^, the third term is still of order 
1/L whereas the second term becomes of order 1. The saddle-point equation to leading 
order in 1/L becomes 

^ = -Trp{z) + ... (8.98) 
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We immediately obtain 



Normalization gives 



We obtain the distribution 
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3 ,^2 2 
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(L^-z^). (8.101) 



We imposed the condition 



^2 o 2 - .37r,l 2 



L'' = 2TTg^a^L={—)3gs. (8.102) 



By comparing with the numerical results we will find that the hermitian prediction is 
more accurate than this antihermitian calculation. 
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9 Summary and Future Directions 

In this article we employed supersymmetry, coliomological deformation, localization 
and the saddle-point method as well as the Monte Carlo method to study nonperturba- 
tively Yang-Mills matrix models in L) = 4 with mass terms. We can summarize the main 
results, findings and conjectures of this work as follows: 

• By imposing the requirement of supersymmetry and 6*0(3) covariance we have shown 
that there exists a single mass deformed Yang-Mills quantum mechanics in D = 4 
which preserves all four real supersymmetries of the original theory although in a 
deformed form. This is the 4 dimensional analogue of the 10 dimensional BMN 
model. Full reduction yields a unique mass deformed D = 4 Yang-Mills matrix 
model. This latter 4 dimensional model is the analogue of the 10 dimensional IKKT 
model. 

• By using cohomological deformation of supersymmetry we constructed a one-parameter 
(Co) family of cohomologically deformed D = A Yang-Mills matrix models which pre- 
serve two supercharges. The mass deformed model is one limit (Co — > 0) of this 
one-parameter family of cohomologically deformed Yang-Mills models. 

• We studied the models with the values /3 = and /5 = 2/9 where j3 is the mass 
parameter of the bosonic matrices Xa- The second model is special in the sense that 
classically the configurations Xa ~ La , ^4 = is degenerate with the configuration 
X„ = 0,^4 = 0. 

• The Monte Carlo simulation of the bosonic Z) = 4 Yang-Mills matrix model with 
mass terms shows the existence of an exotic first/second order transition from a 
phase with a well defined background geometry given by the famous fuzzy sphere 
to a phase with commuting matrices with no geometry in the sense of Connes. The 
transition looks first order due to the jump in the action whereas it looks second 
order due to the divergent peak in the specific heat. 

• The fuzzy sphere is less stable as we increase the mass term of the bosonic matrices 
Xa, i.e. as we increase /3. For j3 = 2/9 we find the critical value a* = 4.9 whereas 
for /3 = we find the critical value (3 = 2.55. 

• The measured critical line in the plane a — /3 agrees well with the theoretical pre- 
diction coming from the effective potential calculation. 

• The order parameter of the transition is given by the inverse radius of the sphere 
defined by 1/r = TrXa/{cP'C2)- The radius is equal to 1/0^ (where (f) is the classical 
configuration) in the fuzzy sphere phase. At the transition point the sphere expands 
abruptly to infinite size. Then as we decrease the inverse temperature (the inverse 
gauge coupling constant) a, the size of the sphere shrinks fast to 0, i.e. the sphere 
evaporates. 

• The fermion determinant is positive definite for all gauge configurations in D = 
4. We have conjectured that the path integral is convergent as long as the scalar 
curvature (the mass of the fermionic matrices) is zero. 
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We have simulated the two models (3 = and /3 = 2/9 with dynamical fermions. The 
model with /3 = has two supercharges while the model /3 = 2/9 has a softly broken 
supersymmetry since in this case we needed to set by hand the scalar curvature to 
zero in order to regularize the path integral. 

Thus /3 = is amongst the very few models (which we known of) with exact super- 
symmetry which can be probed and accessed with the Monte Carlo method. 

The fuzzy sphere is stable for the supersymmetric D = 4 Yang-Mills matrix model 
with mass terms in the sense that the bosonic phase transition is turned into a very 
slow crossover transition. The transition point a is found to scale to zero with N. 
There is no jump in the action nor a peak in the specific heat. 

The fuzzy sphere is stable also in the sense that the radius is equal l/(p'^ over a much 
larger region then it starts to decrease slowly as we decrease the inverse temperature 
a until it reaches at a = 0. We claim that the value where the radius starts 
decreasing becomes smaller as we increase N. 

The model at d = can never sustain the geometry of the fuzzy sphere since it is the 
non deformed model so in some sense the transition to commuting matrices always 
occurs and in the limit N — > oo it will occur at d* — > 0. 

We have spent a lot of time in trying to determine the eigenvalues distributions of the 
matrices X^ in both the bosonic and supersymmetric theories. A universal behavior 
seems to emerge with many subtleties. These can be summarized as follows: 

— In the fuzzy sphere the matrices Xa are given by the SU{2) irreducible rep- 
resentations La- For example diagonalizing the matrix X3 gives N eigen- 
values between (N — l)/2 and —{N — l)/2 with a step equal 1, viz m = 

(N - l)/2, {N - 3)/2, ..., -(iV - 3)/2, -{N - l)/2. 

— In the matrix phase the matrices X^ become commuting. More explicitly the 
eigenvalues distribution of any of the matrices X^ in the matrix phase is given 
by the non-polynomial law 

/'4(x) = ^(i?2-x2)i. (9.1) 

This can be obtained from the conjecture that the joint probability distribution 
of the four matrices X^ is uniform inside a solid ball with radius R. 

— In the matrix phase the eigenvalues distribution of any of the Xa, say X3, is 
given by the above non-polynomial law with a radius R independent of a and 
N. 

— This is also confirmed by computing the radius in this distribution and compar- 
ing to the Monte Carlo data. 

— A very precise measurement of the transition point can be made by observing 
the point at which the eigenvalues distribution of X3 undergoes the transition 
from the TV-cut distribution to the above non-polynomial law. 

— The eigenvalues distribution of X4 is always given by the above non-polynomial 
law, i.e. for all values of a, with a radius R which depends on a and A*". 
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— Another signal that the matrix phase is fuUy reached is when the eigenvalues 
distribution of X/^ coincides with that of X3. From this point downward the 
eigenvalues distribution of X4 ceases to depend on a and N . 

— Monte Carlo measurements seems to indicate that i? = 1.8 for bosonic models 
and R = 2.8 for super symmetric models. The distribution becomes wider in the 
supersymmetric case. 

— We have also observed that the eigenvalues of the normal scalar field X^ — C2 in 
the fuzzy sphere are also distributed according to the above non-polynomial law. 
This led us to the conjecture that the eigenvalues of the gauge field on the back- 
ground geometry are also distributed according to the above non-polynomial 
law. Recall that the normal scalar field is the normal component of the gauge 
field to the background geometry which is the sphere here. 

• In the Z) = 3 Yang-Mills matrix model with mass terms the eigenvalues distribution 
becomes polyonomial (parabolic) given by 

p,{x) = ^^{R^ - x\ (9.2) 

It was difficult for us in this article to differentiate with certainty between the two 
distributions p^ and pa in the three dimensional setting. 

• We have also attempted to compute the above eigenvalues distributions analytically. 
Using localization techniques we were able to find a special set of parameters for 
which the L) = 4 Yang-Mills matrix model with mass terms can be reduced to the 
three dimensional Chern-Simons (CS) matrix model. The saddle-point method leads 
then immediately to the eigenvalues distributions p^. We believe that our theoretical 
prediction for the value of R is reasonable compared to the Monte Carlo value. 

• We have also made a preliminary comparison between the dependence of i? on a in 
the hermitian and antihermitian CS matrix models. The hermitian case seems more 
appropriate for the description of the eigenvalues of X-j, whereas the antihermitian 
case may be relevant to the description of the eigenvalues of X4. 

• Finally, we conjecture that the transition from a background geometry to the phase 
of commuting matrices is associated with spontaneous supersymmetry breaking. In- 
deed mass deformed supersymmetry preserves the fuzzy sphere configuration but 
not diagonal matrices. 

Among the future directions that can be considered we will simply mention the fol- 
lowing four points: 

• Higher precision Monte Carlo simulations of the models studied in this article is the 
first obvious direction for future investigation. The most urgent question (in our 
view) is the precise determination of the behavior of the eigenvalues distributions in 
D = A and Z) = 3. An analytical derivation of pa and especially /94 is an outstanding 
problem. 

• Finding matrix models with emergent 4 dimensional background geometry is also an 
outstanding problem. 
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• Models for emergent time, and to a lesser extent emergent gravity, and as a conse- 
quence emergent cosmology are very rare. 

• Monte Carlo simulation of supersymmetry based on matrix models seems to be a 
very promising goal. 

A Algorithms and Simulations 

We will use the Metropolis algorithm to compute the probability distribution 

dViX^) = ^5 (TrX^) detP e"^^. (A.l) 

In problems involving fermions and/or matrices, the Hybrid Monte Carlo algorithm would 
have been a better choice as it involves global updates as opposed to local updates involved 
in the Metropolis and as a consequence it is more efficient for non local theories. The 
disadvantage in using the Hybrid Monte Carlo algorithm is the fact that the molecular 
dynamics part of this algorithm contains two parameters (the time step size and the 
number of molecular dynamics steps) which require optimization while the Metropolis 
algorithm is free from extraneous parameters. Since we will be reporting data with the 
fermion determinant for small values of A^ up to A^ = 10 we prefer here the Metropolis 
algorithm over the Hybrid Monte Carlo algorithm. For the bosonic theory the Metropolis 
algorithm is quite efficient and we can reach up to A = 100 as in |55y56j. The data with 
the fermion determinant for large values of N using the Hybrid Monte Carlo algorithm 
will be reported elsewhere [203 ]. 

We have used two random numbers in these simulations: RAN2 and RANLUX. We 
have found no discernible or measurable differences between these two generators. For 
consistency, the data reported here is all with RAN2. The errors were estimated using 
the jackknife method which is supposed to take into account and suppress statistical 
correlations. The determinant of the Dirac operator and the eigenvalues are computed 
using LAPACK. 

A simulation consists of Tth + Tmc sweeps where T^h is the number of thermalization 
sweeps and Tmc is the number of Monte Carlo sweeps. For all simulations with the 
fermion determinant we choose T^h = ?mc = 2^'^. For the bosonic simulations we choose 
Tth = T'mc = 2^^ — 2^*^. We have checked that there is no difference between starting 
from a cold start and a hot start in bosonic simulations with /3 = 0. The same holds for 
bosonic simulations with /3 7^ although thermalization becomes somehow difficult. For 
the super symmetric models we always start from the matrices A4 = 0, Xa = (pLa since 
thermalization is very costly in this case. 

In a Monte Carlo sweep we update every single entry of the four matrices X^ once and 
thus a Monte Carlo sweep contains 2(A^ + N) Metropolis steps. This consists one unit 
of time or a step of the Monte Carlo evolution. 

In the actual code there was no symmetry between the matrix X4 and the matrices 
Xa in the sense that we chose different intervals from which we pick the variations of X4 
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and Xa- These intervals are adjusted at each Monte Carlo unit of time (sweep) so that 
the acceptance rates for X4 and Xa are fixed separately between 25 and 30 per cent. 

In actual numerical simulations it is of vital importance to remove the free degrees of 
freedom TrX^, TrOa and Tr9^ in order to achieve a very efficient code. After updating 
the diagonal matrix element i of a given matrix X^ as: {X^)ii — > (^^)m = {X^j)ii + d+d* 
the trace changes as: TrX^ = — > TrX^ = d + d* . Thus to restore a traceless matrix 
X^ we subtract from it the matrix —{d + d*)'\.N /N and as a consequence we change the 
diagonal entries of X^ as 






{X^,)ii = {X^)ii + d + d* 



(X, 



t^!33 



(^, 



M^JJ 



d + d* 

N 



d + d* 

N 



(A.2) 



We found that this is the only prescription of imposing tracelessness that works in the 
matrix phase. The effect of the tracelessness condition on the determinant of the Dirac 
operator is discussed in the next appendix. 



B The Dirac Operator in Simulations 



The aim is to simulate the action 

1 



S 



Sb - J^Tre+ ( i[Xi, ..] + aa[Xa, ..] + a? 



(B.l) 



1 



2ia 



T [^M ' -^<^] "I ~ ^abc^aXbXc 



2t^„ v2 



+ N^a'^TrX^ 



(B.2) 



Sb = NTr 
We work with the effective action 

5eff = 5b -Trad log P. (B.3) 

The Dirac operator is given by 

V = iXi- iX^ + aaXa - cjaX^ + al (B.4) 

It is possible to rewrite the Dirac action in the following form (with X34 = X^ + 1X4) 



TrO'^Ve = Tr 



Tr 



e+{X3i + ai)9i + e+x^e2 + e+x+Oi + e+{-x+^ + ae>2 



X^iO^ Oi + X^O^ 62 + X^62 01 — ^34^2 ^2 



(B.5) 



Following [MUlT], we expand the iVxiV matrices 6*1, 6*2 and 6*5^, 6*^ as 6li = Ea=i ^i^^"^ and 
9+ = J2T=iiGtT{T'^)'^ where the iV x TV matrices T^ are defined by {T^)ij = SuJ 



]3A 
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with iji and j^ are such that A = N{iji — 1) + ja whereas the matrices (T )+ are defined 
by {T )f- = ^iij^Sjjj^ with Ia and ja such that A = N{ja — I) + ia ■ Then we find that 

Tre+V9 = xlMiiXi + X1-M12X2 + X2-^2iX2 + X2-^22X2. (B.6) 

The matrices M-ij are N"^ x N'^ defined by 

{Mu)"^^ = rr(T^)+(X34 + ai)T^ - TrX^^{T^)+T^ 

{M2i)^^ = Tr{T^)+X+T^ - TrX+{T^)+T^ 

{M22)^'' = Tr{T^)+{-X+^ + ai)T'' + TrX+^{T^)+T'' . (B.7) 

We remark that 

Tr{T^)+XT^ - TrX{T^)+T^ = X,,,J,^,, - X,^,J,,,^ 

A = N{JA -l) + iA , B = N{JB - 1) + is. (B.8) 

The A^^— dimensional vectors xi, X2 and xt^ xt ^^^ defined by {xi)A = (^f and {xt)A = 
(of)*- Thus in terms of the 2 A^^— dimensional vectors x and x'^ the above Dirac term 
becomes 

TrO+Ve = x^Mx- (B.9) 

We observe that the trace parts of the matrices X^ drop from the partition function since 
they correspond to Gaussian distributions. Thus the measure should read J dX^5{TrX^) 
instead of simply J dX^. Similarly we observe that if we write 6 = 9q + rjX then the trace 
part 77 will decouple from the rest, viz 



Tre+{i[Xi, ..] + aa[Xa, ..] + a~^\e = Tre+ [i[Xi, ..] + a,[X,, ..] + a~ne^ + aeV^(B.lO) 

Hence the constant fermion modes rja can also be integrated out from the partition function 
and thus we should consider the measure J dOdO^ 5{Tr9a)5{Tr9'^) instead of J dOdO^ . We 
are thus led to consider the partition function 

Z = f dXa 5{TrXa) detP e^^^. (B.ll) 

The determinant is given by 

detP = f dede+6{Tre^)6(Tre+)e^^''^^^^ 



INa ' 



/ dxdx^Sl Y.{Xa)A6^,J^Si Y.(xi)AS^,,^eltx^^X 
•' ^ A=l ' ^ A=\ ' 



I 



I I \ . J I 



dx'dx'^eT^^"-^^^^ . (B.12) 

82 



The vectors Xi-, X2 ^^^ (-^^ ~ 1)— dimensional. The matrix A4. is 2(iV^ — 1) x 2{N'^ — 1) 
dimensional and it is given by 



. / ^/ .' ^' 



We remark that 

M^p^" = ai6^^. (B.14) 

Thus we must have 

detP = detA^'. (B.15) 

We prefer to write D = Dq + a^. The massless Dirac operator Dq will lead to the 
2{N'^ — 1) X 2{N'^ — 1) dimensional matrix Mq defined by 



/ / / / 



(A^o,ii)^-^ = Tr(r^)+X34T^ - TrX34(T^)+r^ 
(A^o,i2)^-^ = Tr{T'^)+X^T^ - TrX^{T^)+T^ 

{Mo,22)^^ = -Tr{T^)+X+^T^ + TrX+^{T^)+T^ . (B.17) 

The 2N'^ X 2A^2 dimensional identity matrix afl2Ar2 will lead to the 2{N'^ - 1) x 2(iV2 _ i^ 
dimensional matrix (x^i^2(N^-i) + t)- This can be calculated as follows. We have 

aiTre+e = a^Y,(^aTOa- (B.18) 

yl=l 

We use ef = - ZaI'i^ ^a^iAjA where A={iA- 1)N + ja- We get 

Af2-lAr2_i 

A=l B=l ^ ^ 

Af2-lAr2_i 

= «^~E Y.^^^y[^■^A^B5^AJB+5^A3A^^B3B)&t■ (B.19) 

A=l B=l ^ ^ 

The first term corresponds precisely to the 2(iV^ — 1) x 2{N'^ — 1) dimensional identity 
matrix whereas the second term defines the matrix 7. Hence we can write 

detP = det(7Wo + a| + a|7). (B.20) 
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C Calculation of the Mass Deformation 

The mass deformed Lagrangian density is 

/:/. = A + 4^1 + ^C2. (c.i) 

3 9 

L.\ — L,^ + •'-mycrs- (^•■^j 



C^ = TrijUaU + ^i?ii7°[y,7'l + ctS't'^V'- 



(C.3) 



^mycrs = ieCijkTrXiXjXk. (C.4) 



C2=Tr(-^SabXaXh\ (C.5) 

We will suppose mass deformed supersymmetric transformations, with a time depen- 
dent parameter e = e{t) which satisfies d^e = /ulle, given by 

^^Xq = (5oXo 
^fiXi = 5oXi 
Sf,i^ = 6oip + /uAe. (C.6) 

The variation of Cq, the Myers term and the fermionic mass term under the new super- 
symmetric transformations is 

^ ^^„7;„0 / 1 r„,0 „,ii r,.v-n ^ u,i „ Ji I v. \- m „-„0^,J r 



<5^£o = ^TrV'7'^(^2f^ '^']^°^*""4[^''^'][^*'^^']""'^^'[^^'^]"^°^ 



^An e. (C.7) 



^^^(Anyers) = ^Tr^i (^ [f , 7^] [^., ^.ItSs) 6. (C.8) 



^S^iTri^M^) = ^Tri^uf- [^'>,Y]DoXi + '-[Y,^^][X„X,] + 2^aV. (C.9) 

The terms containing the covariant time derivative take the form 

■^Trij-f'^Do ( 7°Xn - A - 2MX] e + -^TrijXDoU.e + 2^Tri>-f'^DoM.Xe. (CIO) 
9'' \ J 9^ 9^ 

We have defined M = -j^M-f° = iaU + ^Hijj^H^j^] - 07^7^73. Note that 11, M and 
M are Clifford algebra valued. Requiring the second and third terms to vanish means that 
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n, M and M must be time independent. Requiring the first term to vanish we obtain the 
constraint 

A = 7°Xn-2MX. (C.ll) 

The remainder 



5^ (cq + ^Tr{[,Milj + iee.jk-^TrXiXjXA 



5(2 \^ 4 



^Tr^i - -[y,7^'][X„X,](37°n + 2M - 4c7i7'7^ + 3^67075) 



^l-f'^AU + 2^iMA-^{f,[M,-f^]}[Xi,XJ] e. (C.12) 



We choose 

n = ^7°M - ^(4c + 3e)7°7SS^ = ^l^iaU + 7°^) " ^(6c + 3e)7°7S^7^- (CIS) 
We get 

6^ (co + ^TriPMiP + ieeijk-^TrXiXjXk) = ^TrtP ( - ^7° AH + 2^MA 

- i{Y,[M,^^]}[X^,X,]y. (C.14) 

We can verify that {7*, [M,7-']} = 0. In other words 



S^, Uo + ^TrijMi; + ieeijk-^TrXiXjXk) = ^Tr^j ( - ^7° AH + 2^Ma') e. 

(C.15) 

This must be identified with 

5jc^ + -^TrijMtlj + iee,jk-^TrXiX,Xi\ = -l^Tri;(si,fxX. (C.16) 

In other words we must have 

SijfXj = 7° An - 2MA. (C.17) 

The final result is 

5jco + -^Tri;Mi; + ieei,k-^TrXajXk-^^Si,TrX,X^ = 0. (C.18) 



We compute 



A = -^j'iXH + 3HX) -^iia + ^jVl'')^. (C.19) 
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We get after some algebra 

y 4 o y 

+ f (-a2 + ^e2 + 3za(c+|)y7V)x (C.20) 

We have H = iHijlf,^^ = -iHijeijk^f^j^j'', §[7*, 7^] = <^ijkl^lblk and 7^75 = ^7^7273. 
Thus we can compute that HXH, H^X and XH^ can only be linear in 7*, proportional to 
the cubic product 7^7^73 or proportional to the identity. We can also compute that HX + 
3HX = AHij^^Xj + AHijeijkXkj^j'^j^. The term linear in 7*^ reads —^j^eHijeijkXk- 
This must vanish. Hence either Hij = or e = 0. But since we are interested in a 
non-zero Myers term we take Hij = 0. We get then 

S,,fX, = y(-a2 + ^e2 + 3za(c + |)y7V)x (C.21) 

The cubic product must also vanish thus we get 

c=-^. (C.22) 



In other words we must have 

S.. = fp2 - 

9 



Sij = (e^ - ^a^)6ij. (C.23) 



D Star Products on S^ And R^ and The Flatten- 



ing Limit 



In this appendix we will show how does the star product on the fuzzy sphere approach 
the star product on the Moyal-Weyl plane in the flattening limit. We will follow |196] . 
The commutation relations on the sphere read 

[Xa,Xb] = —=eabcXc- (D-1) 

VC2 



We define the stereographic projections A and A^ in terms of the operators Xa as follows 

A = ^(xi - ix2)B , A+ = Ib{xi + 1x2) , B = — ^— . (D.2) 

2 2 H — X3 



The commutation relation [xi,X2] = iRxj,l ^foi takes now the simpler form 

{A,A^\ = F(\A\^), F{\A\^) = aB 



\^\Af-—\Af- — 



Q = -, \A\^ = AA^. 



(D.3) 



The constraint Yla -^a — ^ reads in terms of the new variables 






^PB' - (/3 + -)i? + 1 + \A\' = 0, f3 = R + a\A\'. (D.4) 
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This quadratic equation can be solved and one finds the solution 

(D.5) 



B^Bi\A\^) = ^ + \ 






a^ a 



For large N we have 4 = j^[l — a'-^]+0{a'^) and hence ^B = ^(l + |^p) or equivalently 
[A,A+] = 2^(l + |^|2)2 + 0(a2). From the formula |^p = L^{^- L^y^L^ it is easy 
to find the spectrum of the operator F(|^p). This is given by 

F{\A\^)\s,m>= F{Xs,m)\s,m> , s = ^-— . (D.6) 



, C2 — m(m + 1) niN — n) . -, -, t,t 

K,m = -ri= T\9 = / ^ , T^ = ^n-i , rn = -s,...,+s , n = s + m+l = l,...,iV. 

[y/c^-m-iy {'^ + s-nY 



Now we introduce ordinary creation and annihilation operators a and a^ which are defined 
as usual by [a, a^] = 1 with the canonical basis \n > of the number operator M = a^a, 
i.e. M\n >= n\n > , n = 0,1, ... , and which also satisfy a\n >= y/n\n — 1 > and 
a~^\n >= \/n + l|n+l >. Next we embed the A^— dimensional Hilbert space H^^ generated 
by the eigenstates \s,m > into the infinite dimensional Hilbert space generated by the 
eigenstates \n > and then define the maps on H]\f given by 

A = fN{N + l)a , A+ = a+fN{N + 1). (D.8) 

It is an identity easy to check that {N + l)ff^{N + I) = \A\'^ and hence (AA + 1)/^(AA + 
l)|n — 1 >= nff^{n)\n — l > and |^p|s,m >= A„_i|s,?tt- >. In other words we can identify 
the first N states | n > in the infinite dimensional Hilbert space of the harmonic oscillator 
with the states of Hj^j via 

|s,m > -H-jn — 1 = s + m > . (D.9) 



As a consequence we have the result /at (n) = w '^^ which clealry indicates that the above 
map is well define (as it should be) only for states n<N . For example A\Q >= fiiN + 
l)a|0 >= because a|0 >= but also because ^|0 >= ^(xi — ix2)B{\A\'^)\s,—s >= 
B{Xs^-s) 2 /~ ^-\^^ ~^ ^~ ^ • The above map also vanishes identically on \s,s >= |A^— 1 > 
as one might check. The relation between F and Jn is easily found to be given by 

F{Xn) = (1 + n)f^{n + 1) - n/^(n). (D.IO) 

Now we define the corresponding coherent states by 



Af-l „ N-1 



Z? 



1 z'^ 

\z;N >= — , > —;= In > , Mn(x) = > , .^ ^ , ,,,,„ 

[/7v(n)]! = /jv(0)/iv(l).../jv(n - l)/iv(n). (D.ll 
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(D.7) 



By construction these states are normalized and they are such that 

1 z^ 

A\z; N >= z\z; N > -^ lA^ - 1 > . (D.12) 

In the large N limit we can check that M^ix) — >{N-l){l+x)^''^ and v^(A^ - l)![/Ar(A^- 
1)]! — >^J'k{N — 1) and hence A\z;N > — >z\z;N > which means that \z;N > becomes 
exactly an ^— eigenstate. 

As it is the case with standard coherent states the above states \z;N > are not or- 
thonormal since < zi;N\z2;N >= Mn{\zi\'^)~2Mn{\z2\'^)~2Mi\i{ziZ2). Using this result 
as well as the completeness relation J diJ,N{z, z)\z; N >< z; N\ = 1 where dfiNiz, z) is the 
corresponding measure we can deduce the identity 

A# (^■^ [^ ( -, Mn{z)Mn{z) 

M7v(l) = / dfiN{z,z) ^ . (D.13) 

This last equation allows us to determine that the measure dfi{z, z) is given by 

d^iv(z, z) = iMN{\z\'^)XN{\z\'^)dz^dz = 2MN{\z\'^)XN{\z\'^)pdpde. (D.14) 

The function X^v is defined by 

27r 

The solution of this equation was found in |196j and it is given by 

iV- 1 
2nXN{x) = Fi{^ + N,^ + N-N + 1--X) , 7 = \/^ —■ (D-16) 

For large A^ where |zp << N we can find that the behaviour of the measure dfj,j\[(z,z) 
coincides with the ordinary measure on S^, viz 

, , _, A^ — 1 idzAdz 
dfiN{z,z)c^ 27r (l + |z|2)2 - ^^-^^^ 

With the help of this coherent state 1^;; A^ > we can therefore associate to every operator 
O a function On{z-, z) by setting < z; N\0\z; N >= On{z-, z). It is therefore clear that the 
trace of the operator O is mapped to the inetgral of the function Om against the measure 
dfij\[{z, z), i.e 

TrO= dpNiz,z)ON{z,z). (D.18) 

Given now two such operators O and P their product is associated to the star product of 
their corresponding functions, namely 

On * Pn{z, -z)^ < z; N\OP\z; N >= f dp{^, f])OMr], ^) .f'^im^'^lfl -P^^^' ^) 

J Mn{\z\^)Mn{\v\ ) 

ONiv,z) = (< z; N\7]; N >)-'^ < z;L\0\r];L > , Pn^z^t]) = (< ri;N\z]N >)-'^ < t]-N\P\z;N > . 

(D.19) 
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The large A^ limit of this star product is given by the Berezin star product on the sphere 
|199j . namely 



A^ - 1 f idrjAdr] 

'2^ . OT 



On * Pn{z, z) = -^^— I j:^-j-—On{ii, z) 



(1 + zrj){l +riz) "' 
(l + |z|2)(l + |r?|2) 



Piv(2,J?).(D.20) 



Finally we comment on the planar (flattening) limit of the above star product which in 
fact is the central point of our discussion here. In this limit we have X3 = —R where 
the minus sign is due to our definition of the stereographic coordinate B in ()D.2p . The 
stereographic coordiantes B, A and A~^ are scaled in this limit as 

S = - , A = —A , A = xi-ix2 , A+ = — i+ , i+ = XI + ix2 , x, = Xi. (D.21) 

This scaling means in particular that the coordinates z and z must scale as z = z/2R and 
z = z/2R. Since we already know that [A, A'^] = ^ /^ (l + |^P)^ in the large A^ limit we 
can immediately conclude that [A,^+] = 29"^ in this limit or equivalently [xi,X2] = —iO"^- 
Next from the result M^ix) — >{N — 1)(1 + x)^~'^ when N — )-oo we can conlude that in 
this large A^ limit we must have 

Mn{\z\^) — ^A^e^l^l'. (D.22) 

The measure dfj,]\f{z,z) behaves as 

^ _ 

dfiN(z,z) = — -^dzAdz. (D.23) 

Putting all these results together we obtain the Berezin star product on the plane [199j . 
namely 

O * P{z, I) = -^ / dfjAdfjOifi, z) e~^^'~^^'-^'^^P{z, fj). (D.24) 

Lastly it is easy to check that the trace behaves in the limit as follows 

^Tr^i^ I dzAl (D.25) 
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Figure 7: The radius of the bosonic models. 
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Figure 10: The average of the Yang-MiUs and Myers actions of the bosonic models. 
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Figure 11: The average of the total action of the bosonic models. 
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Figure 12: The Schwinger-Dyson identity of the bosonic and supersymnietric models. 
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Figure 16: The specific heat of the supersyminetric models. 
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Figure 18: The average of the Yang-Mihs and Myers actions of the mass deformed model. 



113 



The action of the cohomologically deformed model 



(\ m m 



!;sK^*s«*«^^^^s<>K^^^^ *"* * g * m m * M 



CM 



A 
C/} 

V 




-2 
-4 
-6 
-8 
■10 



« • « 



Bosonic,N=16 

N=13 

N=10 

SUSY,N=10 



— I — I 

X J 

X '. 

H H 



N=8 * 
N=6 
N=4 - 



«* 



0.5 1 1.5 



2 2.5 3 3.5 4 
a 



The action of the mass deformed model 



CM 



A 

CO 

V 



2.2 
2 
1.8 
1.6 
1.4 
1.2 



1 1 

Bosonic,N=16 ^ 




1 


1 




1 








N=13 


X ' 














N=10 


X -: 














SUSY,N=10 


B i 














N=8 


--■ 1 












~ 


N=6 


---G 








I * 






N=4 ' 


- -- -•-- -- 








**i 




_ 


* i • 










^-*.* 


*^f 


J 8 ? . i § ^ ^ 
L 


* 'i S 


^ 


1 i i 5 i 


i » 


* ■ 


i 


B \ 


- 


^¥W?¥¥¥¥ 


1 


plPlPSP^ 




xm 




- 


1 1 


1 




1 




- 







3 
oc 



Figure 19: The average of the total action of the supersymmetric models. 
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Figure 21: The eigenvalues distributions of the mass deformed model. 
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Figure 22: The eigenvalues distributions of the mass deformed model. 
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EVs of the N=10 bosonic models with p =2/9 and a =0.5 
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Figure 23: The eigenvalues distributions of X3 of the = 3,4 bosonic models with /3 = 2/9 in 
the matrix phase. 
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EVs of the N=10 bosonic models with p =2/9 and a =10 
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Figure 24: The eigenvalues distributions of the normal scalar fluctuation X^ — C2 of the D = 3,4: 
bosonic models with /3 = 2/9. 
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